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ABSTRACT 


The  development  and  use  of  the  FORTRAN  program  MARGO 
(Marine  Reduction  Gear  Optimization)  is  described.  MARGO 
performs  design  analysis,  weight  minimization,  and  a 
rudimentary  form  of  noise  minimization  using  the  general 
purpose  optimization  program  called  ADS-1  (Automated  Design 
Synthesis,  Version  1) .  Numerous  subroutines  are  presented 
which  calculate  the  associated  design  variables  for  marine 
reduction  gears.  The  entire  program  is  self-documented  and 
easily  modified  by  the  user.  Examples  are  presented  to 
demonstrate  the  utility  of  the  program. 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  gear  design  process  is  an  iterative  one.  Like  the 
"design  spiral"  approach  used  in  naval  architecture,  the 
gear  designer  must  iterate  about  an  initial  design  point, 
assuming  one  value  to  find  another.  With  each  step,  design 
variables  must  be  compared  to  standards  to  ensure  that 
vital  constraints  are  not  violated.  The  process  is  an 
arduous  one  and  constantly  subject  to  human  error.  New 
technologies,  especially  in  the  heat  treatment  of  steels, 
change  the  initial  design  points  and  establish  a  new  set  of 
manufacturing  limitations . 

Perhaps  the  newest  technology  in  machine  design  is 
CAD/CAM  or  "Computer  Aided  Design"  and  "Computer  Aided 
Manufacturing."  Much  has  been  written  about  the  revolu¬ 
tions  sure  to  follow  on  the  heels  of  high  speed  computers. 
While  it  is  true  that  computers  remove  the  engineer  from 
the  drudgery  of  computation  and  especially  from  the  errors 
of  computation,  the  "art"  of  engineering  is  still  required. 
It  is  unlikely  that  a  computer  program  will  be  written 
soon  which  can  design  a  large  machine  without  a  great  deal 
of  human  interaction.  Therefore,  the  reader  is  cautioned 
that  the  program  described  in  this  thesis  is  not  intended 
to  replace  the  gear  designer.  Some  engineers  are  apt  to 
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look  askance  at  any  topic  of  engineering  described  with  the 
word  "computer."  However,  undergraduate  and  graduate 
engineers  are  taught  that  the  advance  of  technology  is  aided 
by  the  use  of  computers .  While  experienced  engineers  have 
learned  from  many  years  of  practice,  the  new  engineer  lacks 
experience  but  has  a  great  deal  of  enthusiasm  for  applying 
computers  to  the  solution  of  engineering  problems. 

The  two  key  components  for  using  a  computer  to  solve 
engineering  problems  of  design  are  logic  and  organization. 
The  experienced  engineer  uses  a  logical  approach  to  solve 
a  problem.  By  formalizing  the  logical  steps  to  solve  a 
problem,  the  computer  can  be  a  useful  tool. 

One  of  the  computer  tools  available  to  engineers  is 
optimization.  Mathematical  in  origin,  research  in  optimi¬ 
zation  has  resulted  in  several  computer  codes  that  have 
proven  themselves  reliable.  Already  widely  used  in  the 
aerospace  industry,  it  is  time  for  marine  engineers  to  apply 
optimization  techniques . 

B,  THESIS  OBJECTIVES 

The  three  most  important  objectives  of  this  thesis 
are  to: 

1.  Apply  the  academic  principles  learned  to  an  actual 
machinery  design  problem. 

2.  Introduce  optimization  techniques  to  the  Navy's 
marine  machinery  design  team. 


3.  Provide  the  Navy's  marine  machinery  design  personnel 
with  an  initial  program  which  illustrates  the  applica¬ 
tion  of  optimization  techniques . 

To  achieve  these  objectives,  the  design  optimization  of 
marine  reduction  gears  was  chosen  as  a  topic.  The  primary 
outcome  of  the  thesis  is  a  program  entitled  "MARGO"  for 
Marine  Reduction  Gear  Optimization.  MARGO  is  a  master 
FORTRAN  program  which  manages  several  groups  of  subroutines 
to  calculate  various  design  parameters  for  double  reduction 
locked  train  gears  like  those  illustrated  in  Figure  1,  used 
in  the  Spruance  class  destroyer  .  MARGO  also  controls  the 
access  to  the  new  generil  purpose  optimization  program  titled 
"ADS-1"  (Automated  Design  Synthesis — Version  1) . 

ADS-1  was  written  by  Professor  Vanderplaats  of  the 
Naval  Postgraduate  School  in  Monterey,  California.  Pro-  . 
fessor  Vanderplaats  is  also  author  of  the  programs  COPES 
(Control  Program  for  Engineering  Synthesis)  and  C0NM1N 
(Constrained  Minimization)  which  are  widely  used  by  NASA 
and  industry.  ADS-1  is  a  new  general  purpose  optimizer  and 
is  described  in  Chapter  II. 

Because  MARGO  is  written  by  a  graduate  student,  there 
may  be  errors  in  -judgement.  Once  again,  the  reader  is  cau¬ 
tioned  that  sound  engineering  judgement  must  still  be 
applied  to  the  factors  which  influence  the  outcome  of  a 
design.  The  adage  "garbage  in,  garbage  out"  is  an  appro¬ 
priate  reminder  and  MARGO  users  are  urged  to  review  carefully 
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the  parameters  used  and  those  omitted.  Nevertheless, 
MARGO  should  serve  as  an  introduction  to  applying  optimi¬ 
zation  techniques  to  marine  machinery  design. 
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II.  OPTIMIZATION  TECHNIQUES 


A.  PURPOSE 

The  purpose  of  this  chapter  is  to  introduce  the  general 
concepts  of  optimization  and  describe  the  optimization 
program  used  in  this  work.  Further  explanations  of  the 
concepts  and  methods  discussed  in  this  chapter  are  found 
in  the  references  mentioned  herein  and  in  texts  on  the  subject 
of  numerical  optimization. 

B .  GENERAL  CONCEPTS 

The  general  nonlinear  constrained  optimization  problem 
is  expressed  in  standard  form  as  follows: 


Minimize  F(X)  (1) 

Subject  to: 

g j  (X)  <  0  j  =  l,m  (2) 

hk(X)  =0  k  =  1,2.  (3) 

xj  <  Xi  <  X^  i  =  l,n  (4) 


F (X)  is  called  the  objective  function  and  represents 
that  property  (e.g.,  weight  or  noise)  which  is  to  be 
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optimized.  Since  most  optimizers  are  written  to  minimize 
the  objective  function,  the  negative  of  F(X)  is  used  when 
it  is  desired  to  maximize  a  property.  For  instance,  if 
f (X)  =  5X  +  10  is  to  be  maximized,  then  for  the  purposes 
of  optimization  the  function  is  expressed  as: 


F  (X)  =  -  f  (X)  =  -  5X  -  10 


F (X)  may  be  a  linear  or  nonlinear  function  of  the  design 
variables  represented  by  the  vector  X.  Further,  the  objec¬ 
tive  function  may  be  implicit  or  explicit  in  X  but  it  should 
be  continuous  and  have  continuous  first  derivatives. 

Equation  (2)  represents  the  inequality  constraints  which 
are  imposed  on  the  design.  For  example,  if  the  stress  is 
not  to  exceed  100,000  psi,  then  the  constraint  would  be 
formulated  as  follows: 


G 


Stress 

100,000 


1  <  0 


A  value  of  stress  less  than  100,000  psi  will  satisfy  the 
constraint  but  a  value  in  excess  of  100,000  psi  will  result 
in  a  positive  G  and  represents  a  violated  constraint. 

Equation  (3)  represents  the  equality  constraints.  For 
instance,  a  specific  length  over  diameter  ( i/d)  ratio  may 
be  desired.  If  so,  the  expression  l  =  1.5d  would  represent 
an  equality  constraint. 


17 


Equation  (4)  represents  side  constraints.  Side  constraints 
are  the  bounds,  upper  and  lower,  within  which  each  design 
variable  must  remain.  For  instance,  a  pinion  facewidth 
shorter  than  10  inches  or  longer  than  30  inches  may  be 
undesirable.  When  a  constraint  is  equal  to  one  of  its  limits 
(say  pinion  facewidth  *  30  inches) ,  that  constraint  is 
said  to  be  active. 

C.  UNCONSTRAINED  MINIMIZATION 

There  are  two  requirements  for  a  solution  to  an  uncon¬ 
strained  minimization  problem: 

i.  The  gradient  of  the  objective  function  (i.e.,  the 

first  derivative  with  respect  to  the  design  variables) 
must  equal  zero,  and 

ii.  The  Hessian  matrix  must  be  positive  definite. 

This  last  requirement  means  that  the  matrix  of  second 
partial  derivatives  of  the  objective  function  with  respect 
to  the  design  variables  (the  Hessian  matrix)  must  have  all 
positive  eigenvalues.  A  positive  definite  Hessian  matrix 
means  that  a  relative  minimum  exists  but  does  not  guarantee 
a  global  minimum  unless  the  Hessian  matrix  is  positive 
definite  for  all  possible  values  of  X.  Unconstrained 
methods  may  be  classified  by  the  order  of  the  derivatives 
required. 

1.  Non-Gradient  Methods 

Random  Search  is  the  simplest  of  all  optimization 
techniques.  The  search  for  a  minimum  is  literally  a  random 
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search  within  the  defined  design  space.  The  combination  of 
randomly  selected  X  vectors  which  produces  the  minimum 
objective  function  is  selected.  The  term  "X  vector"  is 
used  to  describe  a  vector  in  which  "n"  is  the  number  of 
design  variables  used  in  the  problem.  In  order  to  ensure 
that  an  optimum  has  been  located,  the  Random  Search  Method 
requires  a  very  large  number  of  evaluations .  For  this 
reason,  Random  Search  is  the  least  efficient  optimization 
method.  However,  it  is  useful  in  those  circumstances  where 
computer  memory  is  small  and  lends  itself  to  use  on  hand¬ 
held  calculators.  Because  of  the  large  amount  of  machine 
time  required  to  perform  a  random  search,  many  variations 
have  been  devised  to  establish  a  "search  direction"  to 
improve  efficiency. 

Powell's  method  is  probably  the  most  efficient  zero 
order  method  and  was  developed  in  1964  [Ref.  1] .  Based  on 
the  concept  of  conjugate  directions,  Powell's  method  conducts 
an  initial  search  in  n-orthogonal  directions  and  updates 
the  design  by  the  following  equation: 

Xq  =  Xq_1  +  a*Sq  (4) 

where : 

q  *  the  iteration  number; 
a*  *  the  scalar  step; 

S  *  the  vector  search  direction. 
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The  new  search  direction  is  determined  by  connecting 
the  last  design  point  to  the  original  design  point;  thus 
becoming  the  n+1 ,  or  conjugate,  search  direction.  The  method 
breaks  down  when  no  improvement  is  made  to  the  objective 
function  because  subsequent  search  directions  will  not  be 
conjugate.  Also,  after  a  few  iterations  the  search  direc¬ 
tions  begin  to  be  parallel  to  each  other  due  to  numerical 
imprecision.  Powell  overcame  these  problems  with  a  sophis¬ 
ticated  technique  but  the  simplest  solution  is  to  restart 
the  process  with  unidirectional  searches. 

While  there  are  other  zero-order  methods  described 
in  literature  on  optimization,  the  two  methods  described 
above  give  a  suitable  description  of  zero-order  methods. 

The  next  level  of  improvement  is  to  add  gradient  information 
and  use  first-order  methods. 

2 .  Gradient  Methods 

The  simplest  gradient,  or  first-order  method  is  the 
Fletcher- Reeves  Method  of  Conjugate  Directions  [Ref.  2] . 

The  Fletcher- Reeves  method  is  a  modified  form  of  the  steepest 
descent  method  with  an  improved  rate  of  convergence  to  a 
minimum  solution.  The  search  direction  is  selected  according 
to  the  equation: 

Sq  =  -7F (Xq)  +  BqSq-1  (5) 


where : 
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(6) 


Bq  =  tVf(Xq) I2 

!7F(Xq_1)  j2 

V  =  "del"  the  gradient  operator 
-VF(Xq)  =  the  direction  of  steepest  descent 

=  Sq-^  the  previous  search  direction. 

While  the  Fletcher-Reeves  method  offers  improved 
efficiency,  the  Davidon-Fletcher-Powell  and  the  Broydon- 
Fletcher-Goldf arb-Shanno  variable  metric  methods  for  uncon¬ 
strained  minimization  are  more  sophisticated  and  slightly 
more  efficient.  However,  many  engineering  applications  of 
optimization  will  involve  constraints  and  require  a  differ¬ 
ent  algorithm. 

D.  CONSTRAINED  OPTIMIZATION 

Constrained  optimization  techniques  may  be  divided  into 
two  types,  direct  and  indirect  methods.  Direct  methods  treat 
the  constraints  as  limiting  surfaces  of  the  design  space  and 
search  for  an  optimum  within  the  feasible  design  space. 
Indirect  methods  transform  the  constrained  problem  into  a 
set  of  unconstrained  problems. 

1.  Direct  Methods 

Direct  methods  are  popular  constrained  optimization 
algorithms.  One  well  known  direct  method  is  the  method  of 
feasible  directions  [Ref.  3].  The  method  of  feasible 
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directions  determines  a  feasible  search  direction  and  then 


conducts  a  one- dimensional  search  in  that  direction.  This 
one-dimensional  search  continues  as  much  as  possible  without 
violating  a  constraint.  Constraints  which  are  active  are 
used  to  establish  push-off  factors  to  redirect  the  search 
into  the  feasible  region. 

A  recently  developed  direct  method  is  the  Robust 
Feasible  Directions  Algorithm  developed  by  Vanderplaats  in 
1983.  This  algorithm  combines  features  of  the  Method  of 
Feasible  Directions  and  the  Generalized  Reduced  Gradient 
Method.  A  detailed  description  of  this  algorithm  is  presented 
in  Reference  4 . 

2 .  Indirect  Methods 

A  typical  indirect  method  of  constrained  optimization 
is  the  Augmented  Lagrange  Multiplier  method  which  uses 
penalty  functions  to  contain  or  redirect  the  search  into 
the  feasible  design  space.  Even  if  the  initial  design  is 
infeasible,  this  method  will  still  work.  The  Augmented 
Lagrange  Multiplier  method  is  explained  in  more  detail  in 
Reference  5. 

E.  ADS-1  (AUTOMATED  DESIGN  SYNTHESIS— VERSION  1) 

ADS-1  is  a  powerful  purpose  optimization  FORTRAN  program 
for  the  solution  of  nonlinear  constrained  problems.  Funded 
by  a  research  grant  from  the  National  Aeronautics  and  Space 
Administration  (NASA) ,  ADS-1  combines  many  of  the  best 
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optimization  algorithms  into  one  program  with  a  menu  of 
optimization  techniques  from  which  to  choose.  Considerable 
effort  has  been  spent  developing  and  testing  the  algorithms 
presented  in  ADS-1.  Reference  6  is  a  very  thorough  docu¬ 
mentation  of  the  speed  and  accuracy  of  the  majority  of  the 
algorithms  contained  in  a  preliminary  version  of  ADS-1  and 
research  continues  on  the  development  and  improvement  of 
the  program . 

Since  one  of  the  primary  objectives  of  this  thesis  is 
to  introduce  machinery  design  applications  of  optimization 
techniques  and  because  ADS-1  is  the  optimization  program  used 
by  MARGO,  the  following  detailed  information  is  presented  from 
Reference  7  with  permission  of  the  authors. 

1.  Program  Organization 

One  particular  attribute  of  the  ADS-1  program  is  its 
user  friendliness.  Armed  with  a  fundamental  knowledge  of 
optimization,  the  novice  can  begin  to  use  optimization 
immediately.  As  experience  is  gained  and  trust  in  the 
program  develops,  more  and  more  uses  of  optimization  will 
become  apparent.  While  more  experienced  users  will  be  able 
to  select  and  modify  the  more  sophisticated  routines,  it  is 
not  necessary  to  do  so.  To  achieve  this  flexiblity,  the 
program  divides  the  optimization  task  into  three  basic 
levels : 

i.  Strategy--For  example,  Augmented  Lagrange  Multiplier 
Method; 
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ii.  Optimizer-- For  example,  the  Method  of  Feasible 

Directions; 

iii.  One- Dimensional  Search — For  example,  the  Golden 

Section  Method. 

The  ADS-1  program  is  actually  a  subroutine  which 
calls  subordinate  subroutines  to  perform  the  optimization 
task  as  selected  by  the  user.  The  user  simply  passes  the 
desired  options  along  with  other  information  in  the  prescribed 
form  of  the  ADS  subroutine's  calling  arguments.  The  program's 
logic  is  illustrated  in  Figure  2. 

When  the  user  has  set  up  the  initial  parameters 
(described  below)  and  dimensioned  the  required  arrays,  the 
information  parameter,  INFO  is  set  equal  to  -2.  The  ADS 
subroutine  is  then  called.  The  user  is  then  allowed  to 
modify  any  control  parameters,  and  ADS  is  called  again  for 
optimization.  The  ADS  program  initiates  a  search  for  the 
optimum  by  incrementally  adjusting  the  design  variables. 

After  an  initial  adjustment  is  made,  control  is  returned  to 
the  master  program  in  order  to  conduct  an  evaluation  of  the 
objective  and  constraint  functions  which  are  defined  by  the 
user.  The  program  recalls  ADS  and  continues  this  process 
until  no  further  progress  is  made  in  optimization.  At  this 
point  the  ADS  subroutine  has  found  the  optimum,  INFO  is  set 
equal  to  0,  and  the  optimization  is  complete. 

2 .  User  Instructions 

The  calling  statement  and  a  description  of  the 
calling  arguments  is  presented  below: 
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CALL  ADS  INFO,  ISTRAT,  IOPT,  IONED,  NDV,  NCON,  IPRINT, 
IGRAD,  X,  VLB,  VUB,  OBJ,  G,  IDG,  NAC,  IC,  DF,  A, 

NRA,  NCOLA,  WK,  NRWK,  IWK7  NRIWK 


INFO 

Information  parameter.  On  the  first  call  to  ADS, 

INFO  =  -2  or  0,  depending  on  whether  the  user  wishes 
to  over-ride  default  parameters.  On  subsequent 
calls,  when  control  is  returned  to  the  calling  pro¬ 
gram,  INFO  will  have  a  value  of  0 ,  1,  or  2.  If 
INFO  =  0  on  return  from  ADS,  the  optimization  is 
complete.  If  INFO  =  1,  the  user  must  evaluate  the 
objective  and  any  constraint  functions  and  call  ADS 
again.  If  INFO  =  2  the  user  must  evaluate  the 
gradient  of  the  objective  at  a  specified  set  of 
constraints.  If  the  gradient  calculation  control 
is  set  equal  to  zero,  IGRAD  =  0,  all  required 
gradient  information  will  be  calculated  in  ADS  by 
finite  difference. 

ISTRAT 

Optimization  strategy  to  be  used  (Table  1) . 

IOPT 

Optimizer  to  be  used  (Table  2) . 

IONED 

One  dimensional  search  algorithm  to  be  used  (Table  3) . 
NDV 

Number  of  design  variables  contained  in  vector  X. 

NDV  is  the  same  as  n  in  the  mathematical  problem 
statement  in  Equations  (1)  through  (4). 

NCON 

Number  of  constraint  values  contained  in  vector  G. 

NCON  is  equal  to  m+Jl  in  the  mathematical  problem 
statement  given  in  Equations  (2)  and  (3).  NCON  may 
be  zero. 

IPRINT 

A  four-digit  print  control  allowing  varying  levels 
of  output. 
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IGRAD 


Gradient  calculation  control.  If  IGRAD  =  0  is  input 
to  ADS,  all  gradient  computations  are  done  within  ADS 
by  finite  difference.  If  IGRAD  =  1,  the  user  will 
supply  gradient  information  as  indicated  by  the  value 
of  INFO. 


X 


Vector  containing  the  design  variables.  On  the  first 
call  to  ADS,  this  is  the  user's  initial  estimate  of 
the  design,  and  may  or  may  not  define  a  feasible 
design.  On  return  from  ADS,  it  is  the  design  for  which 
function  and  gradient  values  are  required.  On  the 
final  return  from  ADS  (INFO  =  0) ,  X  contains  the 
optimum  design. 

VLB 

A^ray  containing  lower  bounds  on  the  design  variables, 

X  • 

VUB 

Array  containing  upper  bounds  on  the  design  variables, 
Xu. 

OBJ 

Value  of  the  objective  function  corresponding  to  the 
design  defined  by  X.  OBJ  has  the  same  meaning  as  F(X) 
in  the  mathematical  problem  statement  given  in  Equation 
(1). 


G 


Array  containing  the  NCON  constraint  values  correspond¬ 
ing  to  the  current  design,  X. 

IDG 

Array  containing  identifiers  indicating  the  type  of 
constraints  contained  in  G.  Constraints  are  identified 
as  nonlinear  or  linear,  inequality  or  equality. 

NAC 

Number  of  currently  active  and  violated  constraints. 

NAC  is  defined  by  ADS  and  returned  to  the  user. 
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IC 


Array  identifying  currently  active  and  violated  con¬ 
straints  for  which  gradients  are  required.  IC (I) 
gives  the  number  of  the  constraint  (located  in  array 
G) .  Array  IC  is  defined  in  ADS  and  is  returned  to  the 
user. 


DF 


Array  containing  the  gradient  of  the  objective  with 
respect  to  the  current  values  of  X. 


A 


Array  containing  the  gradients  of  the  NAC  constraints 
identified  in  array  IC.  Specifically,  column  J  of 
Array  A  contains  the  gradient  of  constraint  K,  where 
K  =  IC(J) . 

NRA 

Dimensioned  rows  of  A. 

NCOLA 

Dimensioned  columns  of  A. 


WK 


User  provided  work  array  for  real  variables.  WK  is 
used  to  store  internal  scalar  variables  and  arrays 
used  by  ADS . 

NRWK 

Dimensioned  size  of  WK. 

IWK 

User  provided  work  array  for  integer  variables.  IWK 
is  used  to  store  internal  scalar  variables  and  arrays 
used  by  ADS . 

NRIWK 

Dimensioned  size  of  IWK. 

3 .  Strategy  Options 

Table  1  lists  the  strategies  available.  The  param¬ 
eter  ISTRAT  is  sent  to  the  ADS  program  to  identify  the 
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TABLE  1 


ISTRAT 

0 

1 

2 

3 

4 

5 

6 

7 

8 


ADS  STRATEGY  OPTIONS 


_ Strategy  to  be  used _ 

None.  Go  directly  to  the  optimizer. 

Sequential  unconstrained  minimization  using 
the  quadratic  exterior  penalty  function 
method  [Refs.  8,9]. 

Sequential  unconstrained  minimization  using 
the  linear  extended  interior  penalty 
function  method  [Refs.  10-12]. 

Sequential  unconstrained  minimization  using 
the  quadratic  extended  interior  penalty 
function  method  [Ref.  13] . 

Sequential  unconstrained  minimization  using 
the  cubic  extended  interior  penalty 
function  method  [Ref.  14] . 

Augmented  Lagrange  Multiplier  method 
[Refs.  15-19]. 


Sequential  Linear  Programming  [Refs.  20, 

21]  . 

Method  of  Inscribed  Hyperspheres  (Method 
of  Centers)  [Ref.  22] . 

Sequential  Quadratic  Programming  [Refs.  17, 
23,  24]. 
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strategy  the  user  desires.  The  option  of  ISTRAT  *  0  would 
indicate  that  control  should  transfer  directly  to  the 
optimizer.  This  would  be  the  case,  for  example,  when  using 
the  Method  of  Feasible  Directions  to  solve  constrained 
optimization  problems  because  that  optimizer  works  directly 
with  the  constrained  problem.  On  the  other  hand,  if  the 
constrained  optimization  problem  is  to  be  solved  by  creating 
a  sequence  of  unconstrained  minimizations,  with  penalty 
functions  to  deal  with  the  constraints,  one  of  the  appro¬ 
priate  strategies  would  be  used. 

4 .  Optimizer  Options 

Table  2  lists  the  optimizers  available.  IOPT  is 
the  parameter  used  to  indicate  the  optimizer  desired.  The 
option  of  IOPT  =  0  is  not  normally  used.  This  option  is 
provided  for  program  development  where  it  is  desired  to 
access  one  of  the  one-dimensional  search  algorithms  avail- 
able  in  ADS. 

In  choosing  the  optimizer  (as  well  as  strategy  and 
one-dimensional  search)  it  is  assumed  that  the  user  is 
knowledgeable  enough  to  choose  an  algorithm  consistent  with 
the  problem  at  hand.  For  example,  a  variable  metric  opti¬ 
mizer  would  not  be  used  to  solve  constrained  optimization 
problems  unless  a  strategy  is  also  used  to  create  the 
equivalent  unconstrained  minimization  task  via  some  form 
of  penalty  function. 
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0  None.  Go  directly  to  the  one-dimensional 

search.  This  option  is  used  only  for 
program  development. 

1  Method  of  Feasible  Directions  (MFD)  for 
constrained  minimization  [Refs.  3,  25]. 

2  Fletcher- Reeves  conjugate  direction 
algorithm  for  unconstrained  minimization 
[Ref.  2] . 

3  Robust  Method  of  Feasible  Directions  for 
constrained  minimization  [Ref.  4] . 

4  Davidon-Fletcher- Powell  (DFP)  variable 
metric  method  for  unconstrained  minimi¬ 
zation  [Refs.  26,  27]. 

5  Broydon-Fletcher-Goldfarb-Shanno  ( BFGS ) 
variable  metric  method  for  unconstrained 
minimization  (Refs.  28-31]. 

6  Newton's  Method  for  unconstrained 
minimization. 
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5.  One-Dimensional  Search  Options 

Table  3  lists  the  one-dimensional  search  options 
available  for  unconstrained  and  constrained  problems.  Here 
IONED  identifies  the  algorithm  to  be  used.  Normally  the 
option  of  obtaining  bounds  only  should  not  be  used. 

6 .  Allowable  Combinations  of  Algorithms 

Not  all  combinations  of  Strategy,  Optimizer  and 
One-Dimensional  Search  are  meaningful.  For  example,  con¬ 
strained  one-dimensional  search  is  not  meaningful  when 
minimizing  unconstrained  functions. 

Table  4  identifies  the  combinations  of  algorithms 
which  are  available  in  the  ADS  program.  In  the  table, 
an  X  is  used  to  denote  an  acceptable  combination  of  Strategy, 
Optimizer  and  One-Dimensional  Search,  and  it  is  seen  that 
well  over  100  different  algorithms  are  possible.  An  example 
is  shown  by  the  dashed  line  on  the  table  which  indicates 
that  constrained  optimization  is  to  be  performed  by  the 
Augmented  Lagrange  Multiplier  method  (ISTRAT  =  5)  using  the 
BFGS  optimizer  (IOPT  =  5)  and  polynomial  interpolation  with 
bounds  for  the  one-dimensional  search  (IONED  =  4) . 

Because  of  the  vast  number  of  algorithms  developed 
in  recent  years,  it  is  clear  that  this  list  of  options 
could  be  greatly  expanded.  One  of  the  attributes  of  the 
ADS  program  is  the  emphasis  placed  on  its  future  expansion 
with  a  minimum  of  effort.  This  flexibility  is  achieved 
through  modularity  of  the  program  to  avoid  unnecessary 
duplication  of  common  operations. 
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TABLE  3 


ADS  ONE-DIMENSIONAL  SEARCH  OPTIONS 


IONED  One-Dimensional  Search  to  be  used  [Refs.  8,  32] 

1  Find  brackets  on  the  minimum  of  an  unconstrained 
function. 

2  Find  the  minimum  of  an  unconstrained  function 
using  the  Golden  Section  method. 

3  Find  the  minimum  of  an  unconstrained  function 
using  the  Golden  Section  method  followed  by  a 
cubic  polynomial  interpolation. 

4  Find  the  minimum  of  an  unconstrained  function 
by  first  finding  bounds  and  then  using 
polynomial  interpolation. 

5  Find  the  minimum  of  an  unconstrained  function 
by  polynomial  interpolation/extrapolation 
without  first  finding  bounds  on  the  solution. 

6  Find  brackets  on  the  minimum  of  a  constrained 
function. 

7  Find  the  minimum  of  a  constrained  function  by 
using  the  Golden  Section  method. 

8  Find  the  minimum  of  a  constrained  function 
using  the  Golden  Section  method  followed  by 
cubic  polynomial  interpolation. 

9  Find  the  minimum  of  a  constrained  function  by 
first  finding  bounds  and  then  using  polynomial 
interpolation . 

10  Find  the  minimum  of  a  constrained  function  by 

polynomial  interpolation/extrapolation  without 
first  finding  bounds  on  the  solution. 
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TABLE  4 


ADS  PROGRAM  OPTIONS 


Optimizer 


Strategy 

0 

1 

2 

3 

4 

5 

6 

0 

X 

X 

X 

X 

X 

X 

X 

1 

0 

0 

X 

0 

X 

X 

X 

2 

0 

0 

X 

0 

X 

X 

X 

3 

0 

0 

X 

0 

X 

X 

X 

4 

0 

0 

X 

0 

X 

X 

X 

(  c,  ^ 

V 

_ _  (  Y  \ 

y 

[  J  ) - —  —  —  —  - 

^  0 

' 

••“•A* 

_ 

v  a  ; 

It 

6 

0 

X 

0 

0 

0 

0 

•1 

0 

7 

0 

X 

0 

X 

0 

0 

If 

0 

8 

0 

X 

0 

0 

0 

0 

0 

One-Dimensional 

Search 

It 

•1 

1 

X 

0 

0 

0 

0 

0 

II 

0 

2 

X 

0 

X 

0 

X 

X 

If 

X 

3 

X 

0 

X 

0 

X 

X 

II 

X 

4 

X 

0 

X 

0 

X 

(X) 

X 

5 

X 

0 

X 

0 

X 

X 

X 

6 

X 

0 

0 

0 

0 

0 

0 

7 

X 

X 

0 

X 

0 

0 

0 

8 

X 

X 

0 

X 

0 

0 

0 

9 

X 

X 

0 

X 

0 

0 

0 

10 

X 

X 

0 

X 

0 

0 

0 
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7 .  Design  Example 

The  10-bar  truss  shown  in  Figure  3  was  designed 
using  the  ADS  program  to  demonstrate  a  few  of  the  capabili¬ 
ties  available.  This  example  is  a  favorite  design  example 
for  structural  synthesis  [Refs.  33-35]. 

The  cross-sectional  areas  of  10  members  are  taken 
as  independent  design  variables  with  stress  limits  and 
minimum  gage  constraints  imposed  on  each  member.  The  minimum 
gage  is  0.10  square  inches.  The  stress  limits  are  z  25,000 
psi  for  all  members  except  member  9  which  has  a  stress  limit 
of  t  50,000  psi.  The  specific  weight  of  the  material  is 
0.1  pounds  per  cubic  inch  and  the  total  weight  is  to  be 
minimized. 

It  is  recognized  that  direct  optimization  in  member 
space  is  not  the  best  approach  for  structural  optimization 
problems  such  as  this,  where  approximation  techniques  are 
applicable.  Indeed,  even  without  using  approximation  tech¬ 
niques,  a  better  problem  formulation  would  be  to  treat 
reciprocals  of  the  member  areas  as  the  design  variables. 
However,  the  purpose  here  is  simply  to  demonstrate  the 
optimization  capability  of  ADS,  noting  that  design  efficiency 
can  be  greatly  improved  through  careful  problem  formulation. 

The  results  for  various  combinations  of  Strategy, 
Optimizer  and  One-Dimensional  Search  are  listed  in  Table  5, 
but  for  brevity,  all  possible  combinations  are  not  included. 

The  optimum  weights  achieved  are  an  indication  of 


reliability,  whereas  the  number  of  function  evaluations  and 


TABLE  5 


DESIGN  EXAMPLE  OPTIMIZATION  RESULTS 


ISTRAT 

IOPT 

IONED 

Optimum 

Weight 

Number  of 
Analyses 

Number  of 
Gradients 

0 

1 

7 

1516.8 

305 

39 

0 

1 

8 

1532.3 

321 

37 

0 

1 

9 

1519.7 

120 

30 

0 

1 

10 

1530.1 

119 

31 

0 

3 

8 

1497.8 

489 

8 

0 

3 

9 

1497.3 

114 

6 

I 

2 

4 

1648.8 

114 

27 

1 

4 

2 

1534.2 

384 

37 

I 

4 

5 

1549.4 

109 

33 

I 

5 

4 

1575.2 

159 

35 

2 

2 

4 

1603.4 

216 

46 

2 

4 

5 

1522.7 

133 

41 

2 

5 

2 

1505.2 

528 

51 

2 

5 

3 

1635.7 

456 

40 

3 

2 

3 

1619.4 

428 

39 

3 

4 

4 

1511.3 

211 

51 

3 

5 

2 

1505.2 

528 

51 

3 

5 

5 

1713.8 

93 

29 

4 

2 

3 

1617.2 

489 

45 

4 

5 

2 

1505.2 

528 

51 

4 

5 

4 

1500.8 

209 

52 

4 

4 

5 

1718.8 

70 

22 

5 

2 

3 

1527.3 

518 

48 

5 

4 

4 

1504.0 

210 

47 

5 

5 

4 

1496.3 

235 

54 
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gradient  computations  provides  a  measure  of  relative  effi¬ 
ciency.  From  Reference  35,  the  minimum  known  weight  for 
this  problem  is  1497.61  pounds. 

Caution  should  be  exercised  in  drawing  conclusions 
from  one  example  because,  in  general,  both  efficiency  and 
reliability  of  a  given  algorithm  is  problem  dependent. 

Also,  the  results  given  here  are  preliminary  based  on  a 
March,  1983  version.  Program  development  is  continuing  and 
significant  refinements  have  already  been  made. 

•F.  OPTIMIZATION  SUMMARY 

The  foregoing  material  is  intended  to  be  only  an  intro¬ 
duction  to  optimization.  The  body  of  material  written  on 
optimization  is  vast  and  the  curious  reader  can  consult  any 
of  the  numerous  references  for  further  information.  The 
descriptions  included  above  are  deliberately  brief  but  should 
equip  the  reader  to  become  a  user  of  optimization.  While 
research  continues  to  seek  improvements  in  optimizers,  it  is 
sufficient  for  most  preliminary  designs  to  "get  close 
quickly."  As  an  initial  design  becomes  more  refined,  it  may 
be  possible  to  improve  upon  that  design  by  using  a  more 
sensitive  optimizer.  Thus,  a  quick  rough  estimate  could  then 
be  refined  by  a  more  precise  algorithm.  Whenever  the  design 
goal  of  an  engineer  can  be  mathematically  expressed,  that 
goal  can  be  optimized.  Sound  engineering  judgement  is  still 
required  and  with  practice,  an  engineer  can  more  skillfully 
pose  the  problem  for  optimization. 
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As  a  general  purpose,  state-of-the-art,  user-friendly 
optimization  program,  ADS  is  an  excellent  program  to  use. 
The  examples  presented  in  this  work  will  demonstrate  the 
usefulness  of  optimization  techniques  for  the  design  of 
marine  reduction  gears. 


III.  DEVELOPMENT  OF  MARGO 


A.  INTRODUCTION 

Naval  combatants  are  conceived  in  response  to  national 
defense  policies.  Because  policy  is  subject  to  change  and 
interpretation,  so  are  the  designs  of  warships.  While  such 
changes  generate  considerable  consternation  for  the  designer, 
they  are  a  way  of  life  in  the  democratic  process  of  military 
procurements.  Accordingly,  several  sets  of  military  charac¬ 
teristics  may  be  under  consideration  until  the  Congress 
approves  the  final  form  and  funds  are  appropriated  for  a 
new  ship  class. 

These  military  or  combat  characteristics  of  a  warship 
are  first  defined  as  performance  criteria  such  as  maximum 
speed,  endurance,  and  reliability.  The  choice  of  prime 
movers  may  be  nuclear,  gas  turbine,  steam,  diesel,  or  any 
combination  thereof.  The  make-up  of  the  combat  systems  suit 
is  defined  and  redefined  in  response  to  competitive  programs 
and  interests.  Thus,  the  translation  of  combat  characteris¬ 
tics  into  concrete  design  terms  is  not  an  easy  process.  In 
fact,  the  process  is  an  on-going  one  and  subject  to  a  great 
deal  of  criticism  and  frustration  when  late  changes  result 
in  additional  cost  and  delay.  Unfortunately  for  the  military 
designer,  there  is  very  little  control  over  the  political 
forces  which  generate  design  changes.  Rather  than  chafe  at 


these  inevitable  changes,  the  Navy's  design  team.  Naval  Sea 
Systems  Command  must  cope  with  the  situation,  and  the  only 
reasonable  response  is  to  be  flexible. 

Flexibility  is  best  achieved  by  the  careful  development 
of  reliable  computer  algorithms  which  duplicate  the  logical 
thoughts  of  an  experienced  designer.  A  reliable  computer 
algorithm  could  be  used  to  quickly  respond  to  "what  if" 
questions  and  would  reduce  the  natural  frustration  that 
develops  from  the  posing  of  such  questions. 

With  such  thoughts  in  mind,  MARGO  is  presented  as  a 
modest  illustration  of  the  potential  of  computer  programs 
which  combine  the  mathematics  of  optimization  with  the  logic 
of  machinery  design. 

B.  THE  GEAR  DESIGN  PROCESS 

Because  machinery  design  is  an  iterative  process ,  it  is 
difficult  to  describe  the  beginning.  The  usual  procedure  is 
to  assume  one  value  and  derive  several  others  until  something 
detrimental  is  discovered.  However,  this  first  value  must 
originate  from  an  informed  idea  of  what  the  completed  gear 
set  will  be.  In  order  to  have  an  idea  of  what  the  completed 
gear  set  will  be  like,  specifications  are  written  which 
describe  its  completed  form.  The  drafting  of  specifications 
is  subject  to  the  same  uncertainties  described  above,  but 
sooner  or  later  the  specifications  will  be  decided.  In 
the  following  examples,  it  is  assumed  that  a  double  reduction 
helical  gear  set  is  to  be  designed  for  a  50,000  shaft 
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horsepower  (SHP)  guided  missile  destroyer  such  as  the  DDG-51 
class  presently  being  designed.  The  ship  is  to  be  powered 
by  four  General  Electric  LM-2500  marine  gas  turbines  (MGT's) 
with  a  full  power  prime  mover  speed  of  3600  RPM.  Full  power 
RPM  for  the  main  shaft  is  to  be  160  RPM  for  an  overall 
reduction  ratio  of  22.50.  With  these  specifications,  the 
marine  engineer  must  consider  the  following  items  in  the 
design  synthesis. 

1.  Materials 

The  kind  of  steel  to  be  used  is  among  the  first  and 
most  important  decisions  to  be  made.  Material  selection  is 
normally  the  perogative  of  the  manufacturer  in  the  sense 
that  a  design  bid  is  submitted  with  a  particular  steel  in 
mind.  Considerable  attention  must  be  given  to  material 
selection  because  of  its  impact  on  cost,  ease  of  production, 
and  reliability.  For  the  purposes  of  initial  design,  it  must 
be  determined  what  yield  strength  and  hardness  are  required. 
An  upper  limit  is  arbitrarily  set  in  order  to  contain  the 
design  iterations  within  the  range  that  actual  gears  can  be 
manufactured.  Material  constraints  have  long  been  and  will 
continue  to  be  the  area  most  subject  to  change  in  the  design 
of  machinery.  As  harder  and  stronger  materials  become  avail¬ 
able,  machinery  designs  will  be  developed  which  exploit  the 
latest  materials. 

At  present  the  hardness  of  bull  gears  largely  deter¬ 
mines  the  design  possibilities  of  marine  reduction  gears. 
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In  general,  the  harder  and  stronger  the  steel,  the  smaller 
and  lighter  the  gear  box  will  be  for  a  given  SHP  rating. 

Weight  is  a  critical  item  of  interest  in  warship  construction 
and  minimum  weight  is  foremost  on  the  list  of  priorities  for 
a  warship  power  plant.  Less  weight  means  less  cost, 
especially  in  the  form  of  life  cycle  cost  for  fuel.  In 
addition,  the  military  characteristics  of  the  ship  mean  that 
as  much  weight  as  possible  must  be  conserved  in  order  to  fit 
the  desired  combat  systems  configuration  on  the  ship  and  to 
reserve  space  and  weight  for  future  growth.  The  CG-47  class 
cruiser  is  an  excellent  example.  At  9,700  tons  of  displacement, 
this  ship  was  designed  to  be  the  combat  equivalent  of  the 
Russian  Kirov  class  cruiser  which  displaces  22,000  tons. 

Clearly,  U.S.  warships  are  lighter  in  weight  and  smaller  in 
size  than  their  Soviet  counterparts. 

2 .  Gear  Size 

The  size  of  gears  is  the  area  in  which  the  number  of 
possible  combinations  is  infinite.  Ideally,  a  gear  should 
be  sized  to  just  handle  its  maximum  design  load  and  not 
breakdown  until  the  end  of  its  design  life.  Because  of  the 
interest  in  minimum  weight,  modern  marine  gears  are  of 
welded  construction.  Welded  construction  permits  a  wide 
range  of  sizes  and  upper  limits  for  gears  are  determined  by 
the  physical  dimensions  of  heat  treatment  furnaces.  However, 
the  overall  dimension  of  a  gear  box  is  limited  to  what  will 
fit  through  America's  transportation  network.  Unless  the 
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manufacturing  plant  is  located  close  to  a  seaport  and  trans¬ 
portation  is  to  be  by  barge,  railroad  tunnel  clearances 
dictate  the  maximum  height  of  a  gear  box.  Indeed,  the  mili¬ 
tary's  primary  cargo  plane,  the  C-5A  was  designed  to  conform 
to  this  limit  and  13  feet  2  inches  is  the  upper  limit  on 
horizontal  size  for  any  military  materiel  requiring  trans¬ 
portation  by  air  or  rail. 

The  objective  in  gear  sizing  is  to  build  a  gear  large 
enough  and  strong  enough  to  handle  its  full  power  specifica¬ 
tion  and  yet  not  be  excessively  overbuilt.  The  size  is 
determined  by  the  transmitted  torque,  speed,  and  allowable 
stresses.  The  allowable  stresses  in  turn  depend  on  the 
material  used,  its  heat  treatment,  surface  life  and  design 
life.  To  develop  an  improved  design,  multiple  iterations  of 
the  design  are  required. 

The  allowable  surface  compressive  stresses  and  the 
applied  load  will  determine  the  minimum  gear  pitch  diameter 
and  face  width.  The  allowable  bending  stress  will  determine 
the  minimum  tooth  size,  measured  by  diametral  pitch  or  module 
[Ref.  36] .  Allowable  stresses  are  given  in  AGMA  Standards 
[Refs.  37  and  38] . 

In  high  speed  continuous  duty  gear  applications,  such 
as  marine  reduction  gears  and  turbine  drives,  a  very  high 
number  of  load  cycles  may  be  accumulated  on  the  gear  teeth. 

A  warship  operating  at  standard  speed  (15  knots)  with  a  60% 
operations  tempo  for  30  years  would  accumulate  nearly  10^ 
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cycles.  For  these  high  cycle  applications,  it  is  recommended 
that  the  AGMA  allowable  stresses  be  reduced  by  the  factors 
listed  in  Table  6. 


TABLE  6 

ALLOWABLE  STRESS  FACTORS 
No.  Cycles  Allowable  Stress  Factors 
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10 

1.00 

00 

o 

H 

0.90 

io9 

0.81 

1010 

0.729 

io11 

0.656 

io12 

0.590 

The  next  step  is  to  determine  the  required  gear  sizes 
in  terms  of  pitch  diameter,  center  distance,  and  face  width. 
The  Hertz  equation  for  surface  compressive  stress  is 
arranged  as  follows  for  helical  marine  reduction  gears: 

2 

r0.7  TE  (mg+1)  cos  k ,  -i  1/3 

d  i  I  t - : - H  (7) 

t  s  F  mg  sin  v  cos  $  mpj 

C3 


C  =  j  (mg+l) 


(8) 


where : 
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C  =  Center  distance; 
d  =  Pitch  diameter,  pinion; 

E  =  Modulus  of  elasticity  of  gear  material; 

F  =  Facewidth  of  gear  or  pinion; 
mg  =  Gear  ratio  =  D/d; 
mp  =»  Profile  contact  ratio; 

k,  =  Total  derating  factor  (Assume  *  2  for  initial 
designs) 

=  K-K  K 
0  m  v 

=  Overtorque  factor; 

=  Mounting  factor; 

=  Dynamic  Load  Factor; 

T  =  Torque  input  to  pinion; 

■■p  ~  Helix  angle; 

?  =  Normal  pressure  angle 


K 


K, 


K 


m 
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Many  of  the  gear  design  variables  discussed  above  may 
be  unfamiliar.  The  best  source  of  information  on  the  definition 
and  geometry  of  gear  variables  is  AGMA  Standard  115.01  titled 
"Basic  Gear  Geometry."  Appendix  A  contains  descriptive 
figures  of  some  of  the  common  gear  design  variables  and  their 
symbols . 

3 .  Tooth  Size 

After  determining  the  minimum  pinion  and  gear  diameters, 
the  tooth  size  must  be  selected.  For  high  speed  marine  reduc¬ 
tion  gears,  the  general  preference  is  to  use  as  fine  a  pitch 
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as  possible  within  the  limits  imposed  by  tooth  bending  stress. 
This  preference  is  based  on  the  belief  that  minimized  sliding 
velocities  will  reduce  scoring  and  noise. 

However,  it  is  wise  to  size  the  teeth  such  that  the 
critical  stress  used  to  determine  gear  life  will  be  surface 
compression  and  not  tooth  bending.  Using  this  criteria, 
cumulative  fatigue  is  more  apt  to  result  in  pitting  than 
breakage.  Pitting  damage  is  more  likely  to  be  detected  in 
its  early  stages  from  excessive  noise  or  vibration  or  from 
periodic  visual  inspections  of  the  gears  (Naval  gears  are 
visually  inspected  once  a  quarter) .  Should  a  marine  reduction 
gear  fail  first  by  breakage,  the  consequences  are  likely  to 
include  a  complete  failure  of  the  gear  train  which  is  an 
unacceptable  risk.  For  this  reason,  tooth  bending  stress 
is  conservatively  determined. 

The  sizing  of  gear  tooth  proportions  is  a  whole 
science  unto  itself.  Reference  39  discusses  various  popular 
tooth  forms  and  Appendix  A  contains  a  descriptive  figure  of 
gear  tooth  design  variables.  However,  the  following  equations 
are  a  good  starting  point  but  require  considerable  iteration. 

N  «  15[ro?  1  2-~]  cos  $  (9) 

P  mg 


where: 


=  number  of  pinion  teeth. 
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where : 

d  *  Final  pinion  pitch  diameter. 

Determining  the  addendum  also  requires  multiple 
iterations  of  the  equation: 


a  =  1  +  0.5(m  rr  cos  >?t)2(jp  -  jp)  (12) 

p  p  g 

where: 

m  =  Profile  contact  ratio; 

P 

1.65  for  p  -  20 

1.55  for  •*  »  22.5 

1.45  for  5  =  25 

$  =  Transverse  pressure  angle; 

N  =  Number  of  teeth  in  pinion; 

P 

Ng  =  Number  of  teeth  in  gear. 

A  crucial  factor  in  gear  design  is  shaft  size.  In 
practice,  the  shaft  size  drives  the  gear  size.  This  facet 
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further  compounds  the  iteration  process  and  profoundly  impacts 
the  geometric  arrangement  of  the  gear  box.  Manual  means  of 
iteration  are  difficult  because  the  designer  has  difficulty 
grasping  the  sensitivity  of  one  variable's  impact  upon  another. 
For  this  reason,  iteration  by  computer,  especially  by  optimi¬ 
zation  techniques,  is  the  fastest  and  most  reliable  means 
available  to  find  an  optimum  design.  The  optimizer  calculates 
the  necessary  gradient  information  to  steer  the  iterations 
toward  an  improved  design.  Human  iteration  is  often  frustrated 
by  the  inability  to  determine  a  "search  direction"  which  results 
in  an  improved  design. 

C.  PROGRAM  ORGANIZATION 

1.  Brief  Overview 

MARGO  consists  of  a  main  program  which  utilizes 
numerous  supporting  subroutines.  A  data  file  formatted  in 
accordance  with  the  instructions  contained  in  Appendix  B 
supplies  the  master  program  with  required  information  from  the 
user.  Once  the  user  has  formatted  the  data  file,  the  program 
is  executed  by  entering  the  command  word  "MARGO."  Depending 
upon  the  variables  contained  in  the  data  file,  MARGO  will 
perform  weight  minimization  or  noise  minimization.  A  third 
option  is  design  analysis  which  does  not  require  any  calls  to 
the  optimization  program  ADS.  A  flow  chart  of  the  program 
is  illustrated  in  Figure  4. 

2 .  Master  Program 

The  master  program  manages  the  various  logical 


sequences  as  directed  by  the  user  via  the  data  files.  Initial 


no  further  action  is  required  by  the  user  except  to  pick  up 
the  output  for  subsequent  review.  If  the  user  has  access  to 
a  "hard  wired"  terminal  (i.e.,  connected  directly  to  the  main 
frame  computer) ,  the  program  can  be  executed  and  reviewed  on 
the  terminal.  Execution  time  depends  on  the  computer  system 
and  how  many  users  are  logged  on  the  system.  The  examples 
discussed  in  this  work  were  run  on  the  IBM-3033  with  50-100 
users  logged  on  and  execution  times  ranged  from  less  than  10 
seconds  to  3  or  4  minutes.  Longer  execution  times  depend  on 
the  optimization  strategy  selected  and  the  initial  design's 
starting  point  relative  to  the  minimum. 

While  the  program  can  be  run  sitting  at  a  terminal, 
it  is  difficult  to  comprehend  the  output  unless  the  user  is 
interested  in  only  one  variable,  for  example  total  weight. 
Otherwise,  the  5-10  pages  of  output  should  be  carefully 
reviewed  to  ensure  that  the  data  is  reasonable.  There  are 
some  internal  checks  for  unreasonable  values  which  will  be 
highlighted  in  the  output  but  the  user  is  cautioned  to  care¬ 
fully  input  the  data. 

3 .  Supporting  Subroutines 

The  supporting  subroutines  for  MARGO  are  listed  in 
Table  7.  Subroutines  are  numbered  MRGXXX  and  the  number 
sequence  generally  represents  the  sequence  for  use.  Each 
subroutine  is  self-documented  and  normally  narrow  and  specific 
in  purpose.  There  are  no  global  variables  and  each  subroutine 
is  self-contained.  Calling  arguments  are  clearly  delineated 
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.  TABLE  7 


MARGO  SUPPORTING  SUBROUTINES 


Number 

Purpose 

MRG001 

Gear  tooth  number  combinations 

MRG002 

Transverse  Pressure  Angle  (it) 

MRG003 

Pitches  and  Tooth  Proportions 

MRG004 

Line  of  Action  (Z) 

MRG005 

Ratios  and  Factors 

MRG006 

Loads 

MRG007 

Tooth  Data 

MRGO 0  8 

Radius  Calculations 

MRG009 

Involute  Geometry  Factors 

MRGO 10 

Geometry  Factors 

MRGO 11 

Beam  Stress  and  Stress  Concentration  Factors 

MRGO 12 

Compressive  Stresses  and  Tip  Scoring  Factors 

MRGO 13 

Help  Module  (Reserved  for  development  of 
interactive  features) 

MRGO 14 

Aggregate  Contact  Ratio  (for  noise  minimization) 

MRGO 15 

Weight  Estimates 

with  respect  to  input  and  output  values.  The  variables  used 
within  each  subroutine  are  defined  therein.  Variable  names 
conform  with  the  standard  definitions  established  by  the  AGMA 
[Ref.  40]. 

The  subroutine  system  employed  lends  itself  to  modu¬ 
lar  development  and  further  improvement  especially  in  the  area 
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of  improved  "objective  function"  programs.  For  instance, 
reliable  equations  for  noise  minimization  are  apt  to  be  either 
classified  or  proprietary  in  nature.  By  following  the  examples 
presented  for  weight  and  noise  minimization,  combined  with  the 
information  presented  in  Chapter  II  on  the  use  of  ADS,  users 
should  be  able  to  pose  more  complicated  and  reliable  equations 
for  optimization. 

D.  USER  OPTIONS 

MARGO  has  three  user  options  which  are  selected  in 
accordance  with  the  instructions  contained  in  Appendix  B. 

1 .  Design  Analysis 

Option  I  is  similar  to  the  other  two  options  except 
that  the  optimization  program  is  not  called.  The  values  for 
weight  and  aggregate  contact  ratio  (ACR)  are  calculated  but 
the  basic  design  variables  are  not  adjusted.  This  option  is 
useful  for  checking  the  design  analysis  of  the  manufacturer 
or  of  a  previously  manufactured  gear  for  comparison  purposes. 
Design  analysis  calculations  take  less  than  10  seconds  on  the 
IBM-3033  and  are  expected  to  take  less  than  a  minute  on  the 
VAX  Model  11/780  computer. 

The  weight  estimation  subroutine  (MRGC15 )  is  "cali¬ 
brated"  to  reflect  the  weight  of  a  gear  set  completely  con¬ 
structed  of  welded  gears  such  as  the  gear  set  manufactured  for 
the  FFG-7  class.  Thus,  the  weight  estimates  of  older  designs 
suv.  i  as  the  DDG-2  class  will  be  considerably  lighter  than  their 
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actual  weights  because  welded  construction  was  uncommon  on 
earlier  marine  reduction  gears. 

2 .  Weight  Minimization 

Option  II  performs  weight  minimization  using  the 
weight  estimation  subroutine  mentioned  above.  The  "initial 
design"  is  submitted  via  the  data  file  and  all  possible 
design  parameters  are  eligible  for  optimization  within  pre¬ 
scribed  limits.  The  details  and  procedures  for  use  of  this 
option  are  presented  in  the  MARGO  User's  Manual  contained  in 
Appendix  B. 

The  purpose  of  the  weight  estimation  subroutine  is 
to  provide  a  relative  measure  of  design  improvement.  While 
the  estimated  weight  is  probably  within  5%  of  the  manufac¬ 
tured  weight  for  a  welded  gear  set,  the  weight  estimate 
should  be  used  with  caution  for  any  other  application. 

3 .  Noise  Minimization 

Option  III  uses  the  equation  for  aggregate  contact 
ratio  and  is  a  very  rudimentary  prediction  of  radiated  noise. 
Extensive  research  is  underway  to  accurately  predict  noise 
estimates  and  much  of  this  research  is  classified  when  applied 
to  military  applications.  In  order  to  avoid  working  with 
classified  material,  this  work  deliberately  avoided  more 
serious  consideration  of  noise  minimization.  However,  the 
interest  in  noise  minimization  is  probably  more  important  for 
the  design  of  submarine  reduction  gears  and  the  subject  is 
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addressed  in  order  to  attract  the  attention  of  submarine  reduc¬ 
tion  gear  designers  to  optimization  techniques. 

The  details  and  procedures  for  use  of  this  option 
are  also  contained  in  the  MARGO  User's  Manual,  Appendix  B. 
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IV.  TOOTH  COMBINATIONS 


A.  INTRODUCTION 

To  determine  the  tooth  combinations  for  the  design  of 
a  marine  reduction  gear,  the  gear  designer  needs  two 
specifications : 

1)  reduction  ratio 

2)  reduction  ratio  tolerance 

The  first  task  for  the  designer  is  to  find  a  combination 
of  tooth  numbers  that  yield  the  specified  reduction  ratio 
within  the  allowable  tolerance.  The  equation  for  the  overall 
gear  ratio  of  a  double  reduction  gear  set  is: 


R 


N2  v  N4 
NT  N3 


(13) 


where: 

N1  =  Number  of  teeth  on  the  first  reduction  pinion 

N2  =  Number  of  teeth  on  the  first  reduction  gear; 

N3  =  Number  of  teeth  on  the  second  reduction 

pinion; 

N4  =  Number  of  teeth  on  the  second  reduction 

gear . 

In  the  past,  manual  methods  of  selecting  gear  tooth 
numbers  were  complicated  because  large  reduction  gears  for 
marine  propulsion  required  that  each  tooth  on  a  pinion 
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contact  each  tooth  on  a  gear  before  the  same  teeth  mesh 
again.  Known  as  a  "hunting  tooth"  design,  this  requirement 
was  stipulated  in  order  to  prevent  uneven  wear  patterns  from 
developing  on  the  teeth  while  the  gear  teeth  were  cold  worked 
during  initial  operation.  The  litmus  test  for  determining  if 
a  tooth  combination  is  a  "hunting  tooth”  design  is  to  check  for 
common  factors.  A  mathematical  definition  for  no  common 
factors  is  "conjugate"  and  the  terms  "hunting  tooth  design" 
and  "conjugate  set"  are  used  hereafter  to  mean  the  same 
thing. 

For  instance,  the  combination  of  10  and  20  is  not  a 
"hunting  tooth"  combination  because  the  digits  2,  5  and  10 
are  common  factors  of  the  number  20.  However,  the  combination 
13  and  20  is  a  "hunting  tooth"  combination  because  13  is  a 
prime  number  and  not  a  factor  of  20.  It  is  for  this  reason 
that  many  gear  sets  are  manufactured  using  prime  numbers  for 
either  the  pinion  or  the  gear.  However,  there  are  many  other 
possible  combinations  that  will  satisfy  the  "hunting  tooth" 
criteria,  like  9  and  20.  Neither  number  is  a  prime  number, 
but  the  digits  3  and  9  are  not  common  factors  of  20. 

The  difficulty  arose  when  the  designer  had  to  satisfy 
all  three  criteria  at  once;  reduction  ratio,  reduction  ratio 
tolerance,  and  the  hunting  tooth  criteria.  Since  the  choice 
of  tooth  numbers  is  constrained,  considerable  time  was  required 
to  find  a  conjugate  set  by  manual  methods.  However,  these 
stringent  requirements  are  no  longer  required  for  modern 
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marine  propulsion  gears  because  advanced  materials  negate 
the  necessity  for  a  hunting  tooth  design  altogether. 

In  the  past,  bull  gears  were  manufactured  with  Brinell 
hardness  numbers  around  200.  For  materials  this  soft,  cold 
working  would  occur  when  the  gears  were  first  placed  in  opera¬ 
tion.  In  order  to  insure  that  this  cold  working  was  evenly 
distributed,  the  hunting  tooth  design  criteria  was  applied 
except  that  prime  numbers  greater  than  113  were  undesirable 
due  to  manufacturing  limitations. 

Material  science  and  today's  industrial  capability  now 
allow  bull  gears  to  be  hardened  to  Brinell  300-350.  Cold 
working  does  not  occur  with  gears  of  this  hardness  and  the 
hunting  tooth  combination  is  no  longer  required  to  prevent 
cold  working.  However,  there  may  be  other  reasons  for  the 
continued  use  of  the  hunting  tooth  combination  and  the  subrou¬ 
tine  presented  in  this  work  makes  the  use  of  hunting  tooth 
combinations  very  easy  to  attain. 

For  instance,  considerable  research  is  being  conducted 
in  the  area  of  reduction  gear  noise  quieting.  Besides  the 
industrial  interest  in  noise  reduction  for  machinery  there  are 
.many  military  applications  of  noise  quieting  for  reduction  gears 
especially  in  the  field  of  anti-submarine  warfare.  The  number 
of  gear  teeth  meshing  in  a  reduction  gear  plays  a  significant 
role  in  the  amount  of  noise  generated  by  a  reduction  gear. 

In  order  to  provide  a  convenient  computational  means  of 
selecting  and  changing  tooth  combinations,  a  special  subroutine 
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was  written.  The  subroutine  MRG001  finds  a  conjugate  set 
quickly  and  enables  the  designer  to  consider  different  com¬ 
binations.  A  nonconjugate  set  can  be  found  even  faster  for 
those  designs  where  a  hunting  tooth  combination  is  neither 
required  or  desired.  MRG0Q1  also  allows  a  manufacturer  to 
use  a  gear  previously  designed  for  another  application  by 
finding  three  other  tooth  numbers  that  satisfy  all  three 
criteria. 

B.  THE  EFFECT  OF  TOLERANCE  ON  THE  POPULATION  OF  CONJUGATE 

SETS 

The  effect  of  tolerance  on  the  population  of  conjugate 
sets  for  any  given  reduction  ratio  significantly  impacts  the 
time  required  for  the  computer  to  search  and  locate  a 
satisfactory  set.  Yet,  the  reduction  ratio  tolerance  has  a 
less  significant  effect  upon  the  performance  of  the  gear  set. 
For  instance,  the  proposed  reduction  ratio  and  tolerance  for 
the  Navy's  newest  destroyer  design  (DDG-51)  is  22.50  t  0.01. 
The  tolerance  could  actually  be  larger.  Given  a  prime  mover 
speed  of  3600  RPM  and  the  DDG-51's  reduction  ratio  of  22.50, 
the  propeller  shaft  is  turning  160  RPM  at  full  power.  A 
specified  tolerance  of  22.50  ±  0.01  means  less  than  a  0.1 
RPM  change  at  full  power. 


3600 

T2T5T 


159.9289 


3600 

1277? 


160.0711 
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The  full  power  performance  of  a  propeller  is,  for  all  practi¬ 
cal  purposes,  not  affected  by  a  ±1  or  2  RPM  difference. 

Using  this  criteria,  the  reduction  ratio  tolerance  could 
be  increased  to  t  0.14. 


3600 

159 


22.6415 


3600 

161 


22.3602 


Table  8  illustrates  the  effect  of  tolerance  selection 
on  the  number  of  conjugate  sets  or  "hunting  tooth  designs" 
for  a  specified  reduction  ratio  of  22.50.  As  the  numbers 
in  Table  8  illustrate,  the  population  of  conjugate  sets 
declines  with  smaller  tolerance.  Selecting  an  arbitrarily 
small  tolerance  increases  the  search  time  for  the  computer 
to  find  a  set  of  tooth  numbers.  Variable  pitch  propellers 
make  the  reduction  ratio  selection  even  less  restrictive  and 
unless  a  designer  has  a  special  reason  for  specifying  a 
small  tolerance,  the  largest  possible  tolerance  should  be 
allowed. 

Notice  in  Table  8  that  the  results  for  T  =  0.00001  and 
T  =  0.000001  are  identical.  These  results  are  identical 
because  a  finite  number  of  ratios  are  equal  exactly  to 
22.50,  for  example: 


Nl  =  35 

N2  =  78 

N3  =  52 

N4  -  525 


R 


N2  N4 

nT  nT 


78  525 

IT  T T 


22.5  EXACTLY 
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TABLE  8 


EFFECT  OF  TOLERANCE  ON  POPULATION  OF  CONJUGATE  SETS 

FOR  R  =  22.50 


TOLERANCE 

SETS  CONSIDERED 

SETS  TESTED 

CONJUGATE  SETS 

0.1 

2,684,834,520 

8,221,289 

4,477,834 

0.01 

M 

823,444 

306,656 

0.001 

It 

102,604 

30,930 

0.0001 

II 

43,896 

3,199 

0.00001 

n 

42,290 

1,817 

0.000001 

m 

42,290 

1,817 

NOTES : 

"Sets  Considered"  =  The  total  number  of  different  ratios 

( NUM/DENOM )  formed . 

"Sets  Tested"  =  The  total  number  of  different  ratios 

within  the  specified  tolerance. 

"Conjugate  Sets"  =  The  total  number  of  different  ratios 

which  passed  the  tolerance  test 
AND  the  conjugate  test. 

Minimum  Number  of  Teeth  =  NL  =  35 

Maximum  Number  of  Teeth  =  NM  =  850 

Nl  Ranges  From  Nl  -  25%  NM 

N2  Ranges  From  2  <NL  -  NM 
N3  Ranges  From  NL  -  25%  NM 

N4  Ranges  From  7-NL  -  NM 
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The  population  of  conjugate  sets  is  more  than  8.2  million 
for  a  tolerance  of  T  =  0.1.  This  figure  is  more  than  13 
times  as  large  as  the  population  for  T  =  0.01.  It  follows 
that  the  computer  could  find  a  conjugate  set  for  T  =  0.1 
over  13  times  faster  than  for  a  set  where  the  tolerance  is 
specified  for  T  =  0.01.  The  examples  listed  in  Table  9 
further  illustrate  the  speed  of  computation  versus  tolerance 
for  a  specified  reduction  ratio  of  R  =  22.50: 


TABLE  9 

SPEED  OF  COMPUTATION  VERSUS  TOLERANCE 


CPU  Time 


TOLERANCE 

Nl 

N2 

N3 

N4 

RA 

Seconds 

0.1 

35 

71* 

35 

387 

22.4302 

0.17 

0.01 

35 

79* 

35 

349* 

22.5069 

0.33 

0.001 

35 

88 

39 

349* 

22.4996 

30.83 

0.0001 

35 

78 

52 

525 

22.5000 

126.63 

where: 


RA  =  Actual  Reduction  Ratio 
*  Indicates  Prime  Number 

No  Tooth  Number  Was  Specified 


The  examples  presented  above  were  obtained  on  an  IBM- 
3033  computer.  The  IBM- 30 33  is  a  very  large,  fast  computer 
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and  may  not  be  available  to  the  designer.  However,  the  VAX 
Model  11/780  found  the  same  example  set  for  T  =  0.01  in  less 
than  5  CPU  seconds . 

To  show  that  the  results  for  DDG-51's  reduction  ratio 
(R  =  22.50)  are  not  unique,  Table  10  presents  the  population 
sizes  for  several  other  United  States  Navy  warship  reduction 
ratios  for  T  =  0.01. 

C.  PROGRAM  METHOD  DEVELOPMENT 

Reference  41  describes  the  manual  methods  of  finding 
hunting  tooth  combinations  by  conjugate  fraction  methods  and 
proposes  a  computer  algorithm  for  the  process.  However,  the 
example  presented  in  Reference  41  includes  two  sets  which  have 
a  common  factor  of  2  and  the  actual  FORTRAN  program  which 
was  used  is  not  revealed.  Nevertheless,  the  general  method 
proposed  in  Reference  41  was  used  as  a  starting  point  for  the 
development  of  subroutine  MRG001. 

MRG001  sequentially  forms  two  sets  of  integer  pairs. 

The  first  pair  represents  the  product  of  the  tooth  numbers 
for  the  two  pinions  and  the  second  pair  represents  the 
product  of  the  two  gears.  The  equations  used  are: 

NUM  =  N2  x  N4  (13) 

DENOM  =  Nl  x  N3  (14) 


Yielding: 
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TABLE  10 


POPULATION  OF  CONJUGATE  SETS  FOR  UNITED  STATES  NAVY  WARSHIPS 


SETS  SETS  CONJUGATE 


CLASS 

R 

CONSIDE'lED 

TESTED 

SETS 

CGN-36 

18.72 

2,684,834,520 

1,159,764 

427,137 

DD-963 

21.4864 

842,836 

312,459 

DDG-51 

22.50 

823,444 

306,656 

FFG-7 

20.00 

1,037,245 

381,719 

NOTES  : 


R  =  Reduction  Ratio 


"Sets  Considered"  =  The  total  number  of  different  ratios 

(NUM/DENOM)  formed. 

"Sets  Tested"  =  The  total  number  of  different  ratios 

within  the  specified  tolerance. 


"Conjugate  Sets"  =  The  total  number  of  different  ratios 

which  passed  the  tolerance  test 
AND  the  conjugate  test. 


Tolerance 

=  T 

=  0.01 

for  each  case 

Minimum  Number 

of 

Teeth 

=  NL  = 

35 

Maximum  Number 

of 

Teeth 

=  NM  = 

850 
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R  = 


N2  x  N4 
Nl  *  S3 


NUM 

DENOM 


NUM  and  DENOM  are  formed  by  four  nested  "do  loops"  so  that 
Nl ,  N2,  and  N3  are  held  fixed  while  N4  loops  through  all  of 
its  permissible  values.  Then  N2  is  incremented  to  its  next 
value  and  N4  is  looped  through  its  permissible  values  again. 
Once  N2  has  looped  through  all  of  its  permissible  values,  N3 
is  incremented  and  the  process  repeats  itself  until  Nl  is 
incremented  through  the  full  range  of  its  allowed  values. 

A  simple  example  will  illustrate  the  method.  Let  Nl, 

N2,  N3  and  N4  each  range  from  1  to  3;  then  the  iteration 
values  are  as  follows: 


Iteration  Number 


l 


Iteration  Number  Nl 

18  1 

19  2 

20  2 

•  • 

•  • 

35  3 


N3  N2  N4 

3  3  3  ** 

2  1  1  ** 

2  12 

•  •  » 

•  •  * 

•  •  * 

3  3  3 


At  the  position  marked  with  an  asterisk  the  loop 

restarts  at  N4  =  2  because  the  product  for  1  *  2  is  the  same 
as  2  x  l  (the  commutative  law  of  multiplication)  .  By- 
stepping  the  increment  of  N4,  valuable  computation  time  is 
saved  and  the  algorithm  avoids  duplicate  combinations  with 
tooth  numbers  juxtaposed.  For  the  example  above,  this 
feature  provides  a  56%  improvement  in  efficiency  because 
there  are  only  35  different  ratios  (NUM/DENOM)  while  there 

4 

are  3  =  81  different  combinations  for  these  numbers.  For 

the  same  reasons,  the  positions  marked  with  a  double  asterisk 
"**"  have  a  stepped  increment. 

In  actual  practice,  the  magnitude  of  this  efficiency  is 
even  more  pronounced.  The  data  in  Table  8  considered  over 
2.6  billion  different  ratios,  while  the  number  of  different 

4 

combinations  is  equal  to  (850-35)  or  over  441  billion 
combinations.  Without  this  99%  increase  in  efficiency,  an 
IBM- 3033  computer  would  take  over  36  days  of  CPU  time  to 
process  441  billion  ratios.  Any  procedure  which  constrains 
the  area  of  the  search  will  increase  speed.  Additional 
search  speed  efficiency  is  achieved  by  exploiting  the 
observations  described  below. 
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Table  11  illustrates  the  tooth  combinations  and  reduc¬ 
tion  ratios  of  several  warships  in  the  United  States  Navy. 

A  close  look  at  these  gear  sets  reveals  that  first  reduction 
gears  (N2)  have  from  2  to  4  times  as  many  teeth  as  first 
reduction  pinions  (Nl)  have.  Likewise,  bull  gear  teeth  (N4) 
are  apt  to  be  at  least  7  times  the  number  of  second  reduction 
pinion  teeth.  Table  12  contains  the  tabulated  ratios  of  gear 
tooth  numbers  to  pinion  tooth  numbers.  This  point  is  easily 
visualized  by  remembering  that  the  overall  reduction  gear 
ratio  is  the  product  of  the  first  and  second  reductions. 

In  the  case  of  the  DDG-51's  reduction  ratio  of  22.50,  it  is 
clear  that  numbers  around  3  and  7  will  form  a  product  equal 
to  22.50.  The  user  can  exploit  any  knowledge  of  what  the 
general  range  of  first  and  second  reduction  ratios  ought  to 
be  for  the  gear  set  being  designed.  Since  the  first  reduc¬ 
tion  pinion  is  unlikely  to  have  more  than  4  times  the  mini¬ 
mum  number  of  teeth,  the  Nl  loop  should  be  constrained 
between  the  minimum  number  of  teeth  (NL)  and  4*NL.  In  a 
similar  manner,  the  designer  may  decide  that  the  first  reduc¬ 
tion  gear  will  have  no  more  than  25%  of  the  maximum  number 
of  teeth  allowed  (NM) . 

The  minimum  number  of  pinion  teeth  is  constrained  in 
order  to  prevent  undercutting  and  recommended  minimums  are 
published  by  the  American  Gear  Manufacturers  Association 
(AGMA) .  For  instance,  the  AGMA  recommends  tooth  proportions 
for  fine-pitch  involute  spur  and  helical  gears  and  specifies 


TABLE  11 


UNITED 

STATES 

NAVY 

WARSHIP 

1  REDUCTION 

GEARS 

CLASS 

SHP 

N1 

N2 

N3 

N4 

R 

MANUFACTURER 

DD-963 

40 ,000 

99 

235 

58 

525 

21.4864 

Westinghouse 

FFG-7 

40,000 

46 

119 

53* 

410 

20.0123 

Western  Gear 

CGtf-36 

** 

HP  41* 

94 

43* 

351 

18.7154 

General 

Electric 

LP  43* 

99 

CVN-68 

** 

HP  51 

152 

59* 

517 

26.1163 

General 

Electric 

LP  63 

142 

SSN-681 

*★ 

74 

2S5 

99 

776 

30.1884 

Del  aval 

*  Indicates  Prime  Number 
**  SHP  Classified  or  Unavailable 


TABLE  12 

REDUCTION  GEAR-TO-P INION  RATIOS 


CLASS 

N2/N1 

N4/N3 

CGN-36 

HP 

94/41 

2  .  3 

351/43  = 

8.2 

LP 

99/43 

— 

2  .  3 

CVN-68 

HP 

152/51 

= 

3.0 

517/59  = 

8.8 

LP 

142/63 

= 

2.6 

DD-963 

235/99 

= 

2.4 

525/58  = 

9.1 

FFG-7 

119/46 

= 

2.6 

410/53  = 

7.7 

SSN-681 

285/74 

= 

3.9 

776/99  = 

7.8 

YEAR 

1974 

1976 

1971 

1971 

1972 


the  minimum  number  of  pinion  teeth  as  a  function  of  helix 
angle  and  normal  pressure  angle.  The  recommended  minimum 
number  of  fine-pitch  teeth  ranges  from  6  to  32  teeth,  while 
the  minimum  number  of  coarse-pitch  pinion  teeth  for  marine 
propulsion  gear  sets  is  usually  higher. 

The  designer  should  adjust  the  internal  values  of 
MRG001  according  to  the  gear  design  under  consideration. 

For  the  data  presented  in  this  paper,  the  following  values 
were  used: 


Minimum  Number 

of 

teeth 

= 

NL 

= 

35 

Maximum  Number 

of 

teeth 

= 

NM 

as 

850 

Nl  ranges 

from 

NL 

to 

16% 

NM 

(35 

to 

140) 

N3  ranges 

from 

Nl 

to 

16% 

NM 

(Nl 

to 

140) 

N2  ranges 

from 

2xNL 

to 

50% 

NM 

(70 

to 

426) 

N4  ranges 

from 

7xNL 

to 

100% 

NM 

(245 

to 

850) 

Once  MRG001  has  formed  two  product  pairs,  the  ratio  of 
the  pairs  is  compared  to  determine  which  ratios  fall  within 
the  specified  tolerance.  As  mentioned  above,  the  tolerance 
specification  plays  a  major  role  in  determining  the  popu¬ 
lation  size  of  conjugate  sets.  Sets  within  the  specified 
tolerance  are  then  tested  for  conjugacy  by  using  the  FORTRAN 
math  function  "MOD,"  which  returns  the  remainder  of  one 
number  divided  by  another.  For  example,  X  =  MOD  (10,2)  will 
set  X  =  0  because  the  remainder  of  10/2  is  0,  and  Y  =  MOD 
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(11,2)  will  set  Y  =  1  because  the  remainder  of  11/2  is  1. 
MRG001  loops  through  each  integer  from  2  to  the  number  of 
pinion  teeth.  The  loop  begins  at  2  because  MOD  (10/1)  =  0. 
The  zero  value  indicates  that  1  is  a  factor  of  10.  If  10 
is  the  number  of  teeth  in  a  pinion  and  17  is  the  number  of 
teeth  in  the  corresponding  gear,  MRG001  would  have  decided 
that  1  is  a  factor  of  10  and  also  of  17  and  reject  the  pair. 
However,  1  is  a  common  factor  of  every  number,  but  for  the 
purposes  of  meeting  the  hunting  tooth  criteria,  1  is  an 
acceptable  common  factor.  To  prevent  MRG001  from  rejecting 
all  such  combinations,  the  conjugate  test  loops  begin  at  2. 

When  a  hunting  tooth  design  is  required,  MRG001  performs 
a  conjugate  test  on  each  pair  of  teeth  and  both  sets  must 
pass  the  test  to  be  returned  as  a  conjugate  set.  As  soon 
as  a  common  factor  is  found  between  a  pinion  and  gear  set, 
the  pair  is  rejected  and  the  sequential  search  for  another 
set  is  resumed.  Sets  which  pass  the  conjugacy  test  are 
returned  to  the  master  program. 

In  order  to  avoid  consideration  of  any  gear  or  pinion 
with  a  prime  number  greater  than  113,  a  conditional  test  is 
performed  after  each  loop  is  incremented.  Prime  numbers 
greater  than  113  and  less  than  997  are  incremented  by  an 
additional  1.  As  mentioned  above,  manufacturing  limitations 
vis-a-vis  the  numbers  of  teeth  on  a  hobb  preclude  the 
manufacture  of  a  gear  or  pinion  with  these  numbers. 
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V.  COMPARATIVE  DESIGN  RESULTS 

A.  INTRODUCTION 

The  comparison  of  optimized  designs  with  older  reduction 
gears  must  be  considered  carefully.  Modern  techniques  for 
manufacturing  reduction  gears  place  great  emphasis  upon 
conserved  weight.  Thus,  a  comparison  of  the  DDG-2  class 
guided  missile  destroyer's  reduction  gear  (built  in  1960) 
with  an  optimized  modern  design  would  have  little  meaning. 
Hence,  the  comparisons  included  in  this  work  begin  with  the 
DD-963  class,  which  is  considered  a  turning  point  design. 
Designs  prior  to  DD-963  paid  less  attention  to  weight  and 
designs  after  DD-963  were  considered  weight  critical.  The 
FFG-7  class  frigate  reduction  gear  followed  as  a  design  that 
was  clearly  intended  to  be  light-weight.  Finally,  a  sample 
design  for  the  DDG-51  class  guided  missile  destroyer  is 
presented.  Although  DDG-51* s  reduction  gear  is  still  being 
designed,  the  sample  design  included  herein  will  serve  as 
an  interesting  future  comparison  to  the  designer’s  estimate 
of  weight  and  to  the  gear's  actual  manufactured  weight. 

B.  THE  DD-963  CLASS  DESTROYER 

Table  13  contains  a  summary  of  design  variables  for  the 
DD-963  class  destoryer.  Listed  alongside  the  actual  design 
variables  are  the  optimized  counterparts  corresponding  to 
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TABLE  13 


DD-963  REDUCTION  GEAR  DESIGN  VARIABLES 


Variable 

Actual  Value 

Optimized  Value 

General  Data 

SHP 

40,000 

Same 

Reduction  Ratio 

21.4864 

21.4843 

Weight 

167,500 

159,823 

First  Reduction 

Normal  Diametral  Pitch 

6.0 

4.2309 

Transverse  Diametral  Pitch 

5.437847 

3.8345 

Helix  Angle 

25° 

25° 

Pressure  Angle 

14.5° 

14.5° 

Pitch  Line  Velocity 

279.1  ft/sec 

405.6  ft/sec 

Facewidth 

(+  1.87  gap  =  19.12) 

17.25 

29.875 

Center  Distance 

30.7107 

43.552 

Reduction  Ratio 

2 . 3737 

2.3737 

Number  of  Teeth: 

Pinion 

99 

99 

Gear 

235 

235 

Pitch  Diameter: 

Pinion 

18.2057 

25.8185 

Gear 

43.2156 

61.2863 

Second  Reduction 

Normal  Diametral 

Pitch 

4.0 

6.1265 

Transverse  Diametral  Pitch 

3.625231 

5.5525 

Helix  Angle 

25° 

25° 
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Variable 


TABLE  13  (CONT.) 

Actual  Value 


Optimized  Value 


Pressure  Angle 

14.5 

14.5° 

Pitch  Line  Velocity 

103.3  ft/sec 

166.9  ft/sec 

Face  Width 

32.30 

25.33 

(  2.52  gap  = 

35.02) 

Center  Distance 

80.4087 

53.399 

Reduction  Ratio 

9.0517 

9.0508 

Number  of  Teeth 

:  Pinion 

58 

59 

Gear 

525 

534 

Pitch  Diameter: 

Pinion 

15.9990 

10.6259 

Gear 

144.8184 

96.1730 

each  variable.  The  actual  weight  for  the  DD-963's  reduction 
gear  is  167,500  pounds  and  the  optimized  weight  is  159,323, 
a  savings  of  7677  pounds. 

To  arrive  at  this  optimum,  the  MARGO  data  file  was  first 
calibrated  to  reflect  the  material  properties  of  the 
material  used  in  the  DD-963  reduction  gear.  Since  the  de¬ 
tails  of  such  information  is  proprietary  in  nature  and  not 
in  the  public  domain,  the  values  for  design  bending  stress 
and  design  shear  stress  were  manipulated  until  the  MARGO 
weight  estimate  was  fairly  close  to  the  actual  weight.  In 
this  case,  the  MARGO  weight  estimate  was  167,541  pounds  or 
just  41  pounds  over  the  actual  weight.  The  optimization 
program  was  then  executed  using  material  properties  that  are 
quite  similar  to  the  actual  materials. 

The  savings  in  weight  is  accomplished  by  selecting 
design  variables  which  allow  each  component  to  be  stressed 
to  its  upper  limit.  This  situation  should  not  be  disturbing. 
After  all,  an  upper  limit  is  by  definition  the  maximum 
allowed  value  for  a  variable.  Upper  limits  and  appropriate 
K  factors  should  be  adjusted  if  an  additional  margin  of 
safety  and  reliability  is  desired. 

C.  THE  FFG-7  CLASS  FRIGATE 

The  actual  and  optimized  design  variables  for  the  FFG-7 
class  frigate  are  presented  in  Table  14.  The  actual  weight 
of  the  FFG-7  reduction  gear  is  114,168  pounds  and  the  MARGO 
estimate  is  114,370,  or  202  pounds  over  the  actual  weight. 
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TABLE  14 


FFG-7  REDUCTION  GEAR  DESIGN  VARIABLES 


Variable 

Actual  Value 

Optimized  Value 

General  Data 

SHP 

40,000 

Same 

Reduction  Ratio 

20.0123 

20.0027 

Weight 

114,168 

98,262 

First  Reduction 

Normal  Diametral  Pitch 

4.43 

6.8159 

Transverse  Diametral  Pitch 

3.8309 

5.8941 

Helix  Angle 

30.145° 

30.15° 

Pressure  Angle:  Pinion 

6° 

16.0° 

Gear 

16° 

16.0° 

Pitch  Line  Velocity 

163.1  ft/sec 

122.6  ft/sec 

Facewidth 

19.0 

16.4519 

(+  2.75  gap  =  21.75) 

Center  Distance 

18.623 

13.997 

Reduction  Ratio 

2.58696 

2.5870 

Number  of  Teeth : 

Pinion 

46 

46 

Gear 

119 

119 

Pitch  Diameter: 

Pinion 

12.0077 

7.8044 

Gear 

31.0634 

20.1396 

Second  Reduction 


Normal  Diametral  Pitch 

3.5 

6.1074 

Transverse  Diametral  Pitch 

3.0948 

5.4005 

Helix  Angle 

27.840° 

27.84° 
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TABLE  14  (CONT.) 


Variab le 

Actual  Value 

Optimized  Value 

Pressure  Angle 

14.5° 

14.5° 

Pitch  Line  Velocity 

237.9  ft/sec 

162.9  ft/sec 

Facewidth 

(+  3.25  gap  =  31.25) 

28.00 

19.7825 

Center  Distance 

66.143 

45.274 

Reduction  Ratio 

7.73585 

7.7321 

Number  of  Teeth: 

Pinion 

Gear 

53 

410 

56 

433 

Pitch  Diameter: 

Pinion 

Gear 

17.125 

132.4765 

10.3695 

80.1784 

\ 

The  MARGO  weight  estimate  was  obtained  by  using  the  same 
method  described  above  for  the  DD-963.  The  optimized  weight 
is  98,262  or  a  13.9  percent  improvement  of  the  actual  weight. 

While  the  optimized  weights  for  both  the  DD-963  and 
FFG-7  reduction  gears  are  only  slightly  less  than  the  actual 
designs,  each  of  the  optimized  designs  were  attained  in  less 
than  5  CPU  seconds  on  an  IBM- 30 3 3  computer.  Setting  up  the 
data  files  for  the  MARGO  program  took  less  than  ten  minutes 
each.  It  is  likely  that  the  actual  designs  were  obtained 
after  weeks  of  effort  by  several  engineers.  Thus,  MARGO  can 
be  used  to  quickly  locate  a  design  starting  point  for 
further  analysis  by  hand  or  for  more  detailed  analysis  by 
another  computer  code. 

D.  THE  DDG-51  CLASS  GUIDED  MISSILE  DESTROYER 

Since  tne  DDG-51  class  guided  missile  destroyer  is  in 
the  process  of  being  designed,  a  comparison  of  actual  design 
variables  with  optimized  variables  is  not  possible.  Never¬ 
theless,  Table  15  lists  the  optimized  design  variables  for 
subsequent  comparison. 

E.  NOISE  MINIMIZATION  RESULTS 

Noise  minimization  results  were  found  to  be  the  antithe¬ 
sis  of  weight  minimization  results.  In  general,  light¬ 
weight  gears  are  noisier  than  heavier-weight  gears.  This 
is  logical  because  less  energy  is  required  to  vibrate  small 
masses  than  large  masses.  The  most  likely  ise  of  noise 
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TABLE  15 


DDG-51  REDUCTION  GEAR  OPTIMIZED  DESIGN  VARIABLES 

Variable  _ Optimized  Value _ 

General  Data 


SHP 

50,000 

Reduction  Ratio 

22.50299 

Weight 

115,993 

First  Reduction 

Normal  Diametral  Pitch 

3.7399 

Transverse  Diametral  Pitch 

3.2341 

Helix  Angle 

30.15° 

Pressure  Angle 

15.7904° 

Pitch  Line  Velocity 

170.0  ft/sec 

Facewidth 

20.7567 

Center  Distance 

16.388 

Reduction  Ratio 

2.0286 

Number  of  Teeth: 

Pinion 

35 

Gear 

71 

Pitch  Diameter: 

Pinion 

10.8223 

Gear 

21.9538 

Second  Reduction 

Normal  Diametral 

Pitch 

6.710 

Transverse  Diametral  Pitch  5.9333 

Helix  Angle  27.84° 
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Pressure  Angle 
Pitch  Line  Velocity 
Facewidth 


14.1589° 


Center  Distance 
Reduction  Ratio 
Number  of  Teeth: 

Pitch  Diameter: 


113.8  ft/sec 
17.9899 


43.820 

11.0930 

Pinion 

43 

Gear 

477 

Pinion 

7.2472 

Gear 

80.3931 
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minimization  in  conjunction  with  weight  minimization  is  to 
treat  weight  as  a  constraint.  For  instance,  find  the 
quietest  gear  design  that  weighs  no  more  than  120,000  pounds. 
The  optimization  results  will  probably  be  a  gear  designed 
as  close  as  possible  to  the  maximum  allowed  weight.  An 
alternate  method  would  be  to  impose  a  maximum  noise  index 
and  let  the  design  weight  be  as  small  as  possible  for  the 
noise  limitation.  However,  work  in  this  area  depends  upon 
a  more  reliable  means  of  analysis  and  noise  measure. 
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VI.  CONCLUSIONS 


A.  THE  VALUE  OF  OPTIMIZATION  TECHNIQUES 

The  results  in  Chapter  V  illustrate  the  considerable 
time  to  be  saved  by  using  an  optimization  code.  Once  the 
fundamental  program  is  written  for  the  analysis  of  a  gear 
system,  the  optimizer  allows  considerable  flexibility  in 
experimenting  with  design  boundaries.  Probing  the  limits 
of  design  possibilities  via  optimization  is  faster  than 
manual  methods.  Human  manipulation  of  a  multitude  of  design 
variables  is  very  difficult.  An  improving  avenue  of  design 
possibilities  may  be  abandoned  as  the  result  of  a  math 
error.  Computer  programs  are  the  most  reliable  means  for 
performing  repetitive  computations.  By  necessity,  the 
design  process  requires  multiple  iterations,  while  a  human 
may  tire  of  such  repetition,  a  computer  is  unaffected.  Thus, 
the  two  factors  of  error  and  repetition  make  optimization 
techniques  a  desirable  approach. 

B.  MARGO  APPLICATIONS 

While  MARGO  was  designed  for  large  marine  reduction 
gears  in  the  40,000  SHP  range,  it  is  also  applicable  to 
smaller  gear  sets.  By  selecting  the  appropriate  factors, 
MARGO  could  be  used  to  analyze  and  design  double  reduction 
gear  sets  for  generators  and  for  lighter  duty  marine  pro¬ 
pulsion  drives.  Basically  any  double  reduction  helical 
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gear  could  be  designed  and  with  very  minor  modifications, 
spur  gears  as  well. 

C.  AREAS  FOR  FURTHER  DEVELOPMENT 

The  area  most  important  in  military  reduction  gear 
design  is  noise  reduction.  Much  of  the  mathematics  asso¬ 
ciated  with  noise  involves  complex  numbers.  With  an 
increase  in  the  mathematical  difficulty  of  noise  equations, 
the  chances  for  human  error  dramatically  increase.  As  soon 
as  a  credible  analysis  is  found  to  represent  the  noise 
characteristics  of  a  reduction  gear,  that  analysis  should 
be  structured  as  a  subroutine  compatible  with  MARGO.  By 
following  the  example  procedures  for  the  DDG-51  weight 
minimization  problem,  a  sophisticated  noise  minimization 
subroutine  can  be  formulated.  MARGO  can  then  be  modified 
with  the  appropriate  conditional  branches  and  read  state¬ 
ments  to  call  ADS  for  the  optimization.  It  should  also 
be  noted  that  the  examples  presented  above  were  formulated 
as  unconstrained  minimization  problems.  Constraints  on  noise 
limitations  are  a  logical  area  for  future  expansion. 

A  secondary  area  of  improvement  is  to  expand  the  data 
calculated  by  MARGO.  Tooth  geometry  factors  could  be  calcu¬ 
lated  in  greater  detail.  Subroutines  to  design  shafting, 
bearings,  and  lubrication  requirements  are  also  needed.  A 
library  routine  of  past  designs  would  be  useful  for  comparing 
the  trends  of  various  design  variables. 
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ADS  can  be  used  for  any  optimization  problem.  Its 
flexibility  is  demonstrated  by  MARGO,  but  its  capabilities 
are  far  greater  than  the  modest  example  presented.  It  is 
hoped  that  MARGO  has  served  as  an  introduction  to  its 
capability  and  will  inspire  further  applications. 
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APPENDIX  A 


FIGURES  SHOWING  GEAR  DESIGN  VARIABLES 


a. - — 

\ 

\4V 


=  Total  Force 

=  Axial  Component  or  Thrust  Load 
=  Radial  Component 
=  Transmitted  Load 
=  Normal  Pressure  Angle 
=  Transverse  Pressure  Angle 
=  Helix  Angle 


Fiaurs  E .  Basic  Gear  Elements 


addendum 

•  dedenduro 

•  clearance 


7 

fck  *  working  depth 
i(  •  whole  depth 
.  «  circular  pitch 


rf  *  fillet  radius  oi  basic  rack 
t  .  circular  tooth  thickness  -  I 
cp  *  press  tee  angle 


Figure  6 .  invo 


lute  Gear  Teeth 


APPENDIX  B 


MARGO  USER'S  MANUAL 

1.  INTRODUCTION 

Version  1  of  MARGO  is  written  for  use  by  means  of  a 
data  file.  The  user  edits  a  data  file  named  "MARGO  DATA" 
to  enter  the  desired  initial  design  values.  The  program  is 
executed  with  the  command  word  "MARGO."  The  MARGO  exec 
is  written  and  installed  at  the  Naval  Sea  Systems  Command 
Design  Automation  Center  to  perform  the  utility  commands 
associated  with  compiling  and  running  the  program.  After  the 
MARGO  command  word  is  entered,  the  program  is  executed 
without  any  further  input  from  the  user.  Printed  output 
can  be  picked  up  at  the  user's  designated  pick-up  point  for 
subsequent  review. 

While  Chapter  II  discusses  the  use  of  ADS,  MARGO  Version 
1  does  not  allow  the  user  to  externally  adjust  or  select 
any  of  the  optimization  parameters  described  therein.  The 
optimization  options  used  within  MARGO  are  described  in 
Section  5  of  this  appendix.  Users  who  desire  to  experiment 
with  different  optimization  strategies  should  make  a  copy  of 
the  MARGO  program  under  another  name  and  modify  the  optimi¬ 
zation  segment  to  suit  the  user's  needs.  MARGO  and  its 
supporting  subroutines  are  self-documented  and  modification 
should  be  easy  with  a  printed  copy  of  the  program  to  edit. 
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MARGO  reserves  space  for  an  interactive  option.  How¬ 
ever,  this  option  was  not  developed  in  Version  1  in  view 
of  the  limited  access  to  terminals  directly  connected 
("hard-wired")  to  the  main  computer  at  the  Naval  Sea  Systems 
Command  Design  Automation  Center.  Interactive  operation  on 
a  remote  terminal  with  less  than  a  1200  baud  data  transfer 
rate  is  undesirable.  Initial  versions  of  the  program  were 
run  in  an  interactive  mode  on  an  IBM- 30 33  computer  and  com¬ 
pared  to  the  VAX  Model  11/780.  The  IBM  runs  considerably 
faster  than  the  VAX,  but  even  on  the  VAX  the  analysis  program 
can  be  run  in  less  than  1  CPU  minute. 

2.  HOW  TO  USE  MARGO 

An  initial  data  set  is  provided  for  the  DDG-51.  This 
data  set  should  be  copied  and  stored  under  a  new  name  such 
as  " DDG 5 1  DATA."  The  "MARGO  DATA"  file  can  be  edited  and 
changed  to  the  user's  preference.  Once  the  data  set  is 
fixed,  save  the  file  and  execute  the  program  by  entering  the 
command  word  "MARGO."  The  data  is  echo-printed  in  the  output 
under  both  print  output  options  (Detailed  and  Summary  print 
options).  Thus,  the  key  to  using  Version  1  is  editing  the 
data  file  to  reflect  the  desired  values.  The  data  file 
values  are  explained  below. 

3.  DATA  FILE 

The  data  file  consists  of  30  data  lines  described  below 
and  summarized  in  Table  16.  A  sample  data  file  is  presented 
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in  Figure  8.  Additional  data  lines  may  be  added  by  the 
user  and  inserted  in  Segment  0000,  INITIAL  DATA  INPUT  of 
the  master  program.  Integer  data  are  in  format  110  and  real 
data  are  in  format  F15.5.  This  system  makes  a  sight-check 
of  the  data  file  easy. 

Data  Line  Number  1:  ITCON 

ITCON  =  Iteration  Control.  This  variable  determines  the 
iteration  option  of  the  user  according  to  the  following 
single  integer  code  placed  in  column  10 . 

0  =  Interactive  mode  (reserved  for  subsequent) 

modifications) 

1  =  Analysis  only 

2  =  Optimization  and  Analysis  (of  the  optimum) 

Data  Line  Number  2 :  IOC 

IOC  =  Integer  Optimization  Code.  This  variable 
determines  the  optimization  path  desired  according  to 
the  following  single  integer  code  placed  in  column  10. 

0  =  No  Optimization  (use  this  code  for  analysis 

only) 

1  =  Noise  minimization 

2  =  weight  minimization 

Data  Line  Number  3:  IPC 

IPC  *  Integer  Print  Control.  IPC  determines  the  level 
of  printed  output  desired  according  to  the  following 
single  integer  code  placed  in  column  10. 

1  =  Summary  Output  Only 

2  =  Detailed  Output 

Data  Line  Number  4-7:  Nl,  N2 ,  N3,  N4 

Nl  =  number  of  teeth  on  the  first  reduction  pinion. 

N2  =  number  of  teeth  on  the  first  reduction  gear. 

N3  =  number  of  teeth  on  the  second  reduction  pinion. 

N4  =  number  of  teeth  on  the  second  reduction  (bull) 

gear. 

When  ITCON  =  1,  analysis  is  performed  on  the  designated 
design  for  the  number  of  teeth  specified  and  the 
conjugate  tooth  finder  subroutine  is  not  called.  When 
ITCON  =  2,  the  conjugate  tooth  finder  subroutine  MRG001 
is  called.  Any  single  tooth  number  can  be  specified  or 
a  complete  tooth  set  can  be  specified  and  MRG001  will 
find  the  next  sequential  set.  See  Chapter  IV  for 
additional  information  about  the  use  of  subroutine 
MRG001.  Tooth  numbers  must  be  right  justified  beginning 
in  column  10.  Placing  zeros  in  column  10  will  cause 
MRG001  to  find  a  proper  tooth  combination  for  design. 
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Data  Line  Number  8-9:  AAD(l) ,  AAD ( 2 ) 

AAD(l)  =  the  helix  angle  for  the  first  reduction  and 
AAD (2)  =  the  helix  angle  for  the  second  reduction.  When 
ITCON  =  1,  analysis  is  performed  according  to  the 
specified  helix  angle.  When  ITCON  =  2,  helix  angle  is 
treated  as  a  design  variable  by  ADS  and  will  be 
adjusted  to  its  optimum  value.  Values  are  written  in 
degrees  with  the  decimal  point  in  column  10. 

Data  Line  Number  10-11:  APND(l) ,  APND(2) 

APND(l)  =  the  normal  pressure  angle  for  the  first 
reduction  mesh  and  APND(2)  =  the  normal  pressure  angle 
for  the  second  reduction  mesh.  When  ITCON  =  1,  analysis 
is  performed  according  to  the  specified  pressure 
angles.  When  ITCON  =  2,  normal  pressure  angle  is 
treated  as  design  variable  by  ADS  and  will  be  adjusted 
to  its  optimum  value.  Values  are  written  in  degrees 
with  the  decimal  point  in  column  10. 

Data  Line  Number  12-13:  BL(1) ,  BL(2) 

BL ( 1)  =  the  backlash  for  the  first  reduction  mesh  and 
BL (2)  =  the  backlash  for  the  second  reduction  mesh. 
Normal  values  range  from  0.010-0.020  inches.  Values  are 
written  with  the  decimal  point  in  column  10. 

Data  Line  Number  14:  BUL 

BUL  =  Bearing  Unit  Load  and  is  used  in  the  weight 
estimation  subroutine  MRG015.  A  normal  value  is  300 
psi.  Write  the  value  as  a  real  number  with  the  decimal 
in  column  10. 

Data  Line  Number  15-16:  F(l),  F(2) 

F ( 1)  =  the  facewidth  of  the  first  reduction  mesh  and 
F ( 2)  =  the  facewidth  of  the  second  reduction  mesh.  When 
ITCON  =  1,  analysis  is  performed  using  the  indicated 
values.  When  ITCON  =  2,  ADS  will  adjust  facewidth  as  a 
design  variable  and  determine  its  optimum  value.  Values 
are  written  as  real  numbers  with  the  decimal  point  in 
column  10. 

Data  Line  Number  17:  FLT 

FLT  =  the  locked  train  factor.  Since  MARGO  is  written 
for  double  reduction  gear  sets,  set  FLT  =  2.0  with  the 
decimal  point  in  column  10. 

Data  Line  Number  18:  HP 

HP  =  the  horsepower  to  be  transmitted  to  the  main  shaft 
(SHP) .  Write  this  value  as  a  real  number  with  the 
decimal  point  in  column  10. 

Data  Line  Number  19-20:  P(l),  P(2) 

P(l)  =  the  diametral  pitch  of  the  first  reduction  mesh 
and  P ( 2 )  =  the  diametral  pitch  for  the  second  reduction 
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mesh.  When  ITCON  =  1,  analysis  is  performed  using  the 
indicated  value.  When  ITCON  =  2,  ADS  treats  diametral 
pitch  as  a  design  variable  and  will  adjust  its  value  to 
the  optimum.  Write  diametral  pitch  as  a  real  number 
with  the  decimal  in  column  10. 

Data  Line  Number  21:  PT 

PT  =  propeller  thrust  and  is  used  in  the  weight 
estimation  subroutine  MRGQ15.  A  value  of  500,000  pounds 
is  typical  for  design  purposes.  Write  this  value  as  a 
real  number  with  the  decimal  in  column  10. 

Data  Line  Number  22:  R 

R  =  the  specified  reduction  ratio  when  designing  and 
the  actual  reduction  ratio  when  analyzing. 

Data  Line  Number  23-24:  RF(1),  RF(2) 

RF ( 1)  =  tooth  fillet  radius  for  the  first  reduction 
mesh  and  RF(2)  =  tooth  fillet  radius  r  the  second 
reduction  mesh.  A  normal  value  for  d  oign  purposes  is 
0.02  inches.  Write  these  values  as  real  numbers  with 
the  decimal  in  column  10 . 

Data  Line  Number  25:  RPM 

RPM  =  the  revolutions  per  minute  for  the  prime  mover 
driving  the  reduction  gear.  For  the  LM2500  marine  gas 
turbine,  this  value  is  3600  rpm.  Write  this  value  as  a 
real  number  with  the  decimal  point  in  column  10. 

Data  Line  Number  26-27:  RT(1),  RT(2) 

RT ( 1)  =  the  edge  radius  of  the  generating  rack  for  the 
first  reduction  mesh  and  RT(2)  =  the  edge  radius  of  the 
generating  rack  for  the  second  reduction  mesh.  Write 
these  values  as  real  numbers  with  the  decimal  in  column 
10. 


Data  Line  Number  28:  SDB 

SDB  =  design  bending  stress.  Use  12,500  psi  as  a  normal 
design  point.  Write  this  value  as  a  real  number  with 
the  decimal  point  in  column  10. 

Data  Line  Number  29:  SDS 

SDS  =  design  shear  stress.  Use  7,500  psi  as  a  normal 
design  point.  Write  this  value  as  a  real  number  with 
the  decimal  in  column  10. 

Data  Line  Number  30:  7 

T  =  the  reduction  ratio  tolerance  specified  for  an 
initial  design.  See  Chapter  IV  for  additional 
information  about  specifying  T.  This  value  is  only  used 
when  the  tooth  combination  subroutine  MRG001  is  called. 
Write  this  value  as  a  real  number  with  the  decimal 
point  in  column  10. 
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TABLE  16 


MARGO  DATA 

FILE  LINE  DESCRIPTIONS 

Data  Line 
Number 

MARGO 

Variable 

Remarks 

1 

ITCON 

Iteration  Control 

2 

IOC 

Integer  Optimization  Code 

3 

IPC 

Integer  Print  Control 

4 

N1 

First  Reduction  Pinion  Teeth 

5 

N2 

First  Reduction  Gear  Teeth 

6 

N3 

Second  Reduction  Pinion  Teeth 

7 

N4 

Second  Reduction  Gear  Teeth 

8 

AAD(l) 

First  Reduction  Helix  Angle 

9 

AAD  ( 2 ) 

Second  Reduction  Helix  Angle 

10 

APND(l) 

First  Reduction  Pressure  Angle 

11 

APND ( 2 ) 

Second  Reduction  Pressure  Angle 

12 

BL(1) 

First  Reduction  Backlash 

13 

BL  ( 2 ) 

Second  Reduction  Backlash 

14 

BUL 

Bearing  Unit  Load 

15 

F  ( 1 ) 

First  Reduction  Facewidth 

16 

F(2) 

Second  Reduction  Facewidth 

17 

FLT 

Locked  Train  Factor 

18 

HF 

Horsepower 

19 

P(l) 

First  Reduction  Pitch 

20 

P  (2) 

Second  Reduction  Pitch 

21 

PT 

Propeller  Thrust 

22 

R 

Reduction  Ratio 

23 

RF  ( 1 ) 

First  Reduction  Tooth  Fillet  Radius 

24 

RF  ( 2 ) 

Second  Reduction  Tooth  Fillet  Radius 

25 

RPM 

Prime  Mover  RPM 

26 

RT  ( 1) 

First  Reduction  Edge  Radius  of 
Generating  Rack 

27 

RT  ( 2) 

Second  Reduction  Edge  Radius  of 
Generating  Rack 

28 

SDB 

Design  Bending  Stress 

29 

3DS 

Design  Shear  Stress 

30 

T 

Reduction  Ratio  Tolerance 
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Note:  Integer  data  is  in  format  110  and  real  data  are  in 

format  F15.5.  This  system  makes  a  sight-check  of  the  data 
file  easy  (i.e.,  the  decimal  points  line  up  with  single 
digits)  . 

4.  MARGO  ORGANIZATION 

MARGO  is  self-documented  to  the  extent  that  each 
variable  is  defined  within  each  subroutine  and  all  variables 
are  defined  within  the  master  program.  In  addition,  the 
master  program  is  divided  into  segments  which  are  clearly 
identified. 

The  segment  numbers  and  their  contents  are  defined  as 
follows : 

Segment  0000  Initial  Data  Input 

Segment  1000  Data  File  Echo-Print 

Setment  2000  Interactive  Option  (reserved  space) 

Segment  7000  Optimization 

Segment  7100  Noise  Minimization 

Segment  7200  Weight  Minimization 

Segment  7300  Noise  and  Weight  Minimization 

(reserved  space) 

Segment  8000  Help  (reserved  space) 

Segment  9000  Printer  Output  Options 

5.  MARGO  ADS  PARAMETERS 

MARGO's  ADS  parameters  for  weight  minimization  are 
imbedded  in  the  main  program  in  segments  7000  and  7200. 


These  parameters  must  be  changed  within  the  main  program. 
Later  versions  of  MARGO  can  be  modified  to  permit  external 
changes  to  be  selected  values  at  the  option  of  the  user. 
However,  this  option  was  omitted  in  Version  1  in  order  to 
reduce  the  burden  on  beginning  users.  The  design  variables 
used  in  the  program  for  weight  minimization  are  summarized 
in  Table  17  and  their  ADS  translations  are  listed  in  Table 
18. 


TA3LE  17 

DESIGN  VARIABLES  FOR  WEIGHT  MINIMIZATION 


Design  Element 

MARGO  Variable 

Lower  Bound 

Upper  Bound 

Facewidth 

F  ( 1 ) 

10.0 

24.0 

F  ( 2 ) 

18.0 

36.0 

Pressure  Angle 

APND(l) ,  APND ( 2 ) 

5 . 0 

25.0 

Helix  Angle 

AAD(l) ,  AAD ( 2 ) 

10.0 

35.0 

Pitch 

P(l) ,  P(2) 

2.0 

8.0 

MARGO's  ADS  parameters  for  noise  minimization  are 
imbedded  in  Segments  7000  and  7100  of  the  main  program. 
Because  the  noise  minimization  objective  function  is  not 
considered  to  be  a  reliable  measure  of  noise,  only  the  first 
reduction  mesh  is  analyzed.  This  procedure  is  considered 
an  acceptable  indication  of  relative  noise  measure  since 
the  first  reduction  mesh  operates  at  higher  speeds  than  the 
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TABLE  18 


ADS  PARAMETERS  FOR  WEIGHT  MINIMIZATION 


ADS  Parameter 

MARGO  Value 

NRA 

10 

NCOLA 

5 

NRWK 

1000 

NRIWK 

800 

INFO 

0 

ISTRAT 

0 

IOPT 

5 

IONED 

4 

IPRINT 

2220 

IGRAD 

0 

NDV 

8 

NCON 

0 

X(l) 

F  ( 1) 

X(2) 

F(2) 

X  ( 3) 

APND(l) 

X  ( 4) 

APND(2) 

X(5) 

AAD(l) 

X(  6) 

AAD(2) 

X(7) 

P(l) 

X{8) 

P(2) 

VLB (1) 

10.0 

VLB (2) 

18.0 

VLB (3) 

5.0 

VLB (4) 

5.0 

VLB (5) 

10.0 

VLB (6) 

10.0 

VLB ( 7) 

2.0 

VLB (8) 

2.0 

VUB(l) 

24.0 

VUB (2) 

36.0 
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TABLE  18  (CONT.) 


ADS  Parameter 

VUB (3) 

VUB(4) 

VUB (5) 

VUB (6) 

VUB (7) 

VUB (8) 

OBJ 


MARGO  Value 

25.0 

25.0 

35.0 

35.0 

8.0 

8.0 

TW  (Total  Weight) 


second  mesh.  For  this  reason,  the  first  reduction  mesh 
generates  more  ambient  noise  than  the  second  reduction  mesh. 

This  argument  ignores  the  effects  of  harmonics  which 
carry  farther  in  water  for  lower  frequency  sounds  such  as 
those  generated  by  the  second  reduction  mesh.  However,  a 
detailed  discussion  of  the  noise  characteristics  of  reduc¬ 
tion  gears  with  respect  to  underwater  sound  detection  is 
not  possible  in  this  work  because  of  the  classified  nature 
of  such  subjects.  Thus,  the  noise  measure  minimized  is 
ambient  noise. 

The  noise  minimization  design  variables  are  summarized 
in  Table  19,  and  Table  20  lists  the  ADS  parameters  used. 

TABLE  19 

DESIGN  VARIABLES  FOR  NOISE  MINIMIZATION 


Design  Element 

MARGO 

Variable 

Lower 

Bound 

Upper 

Bound 

Effective  Facewidth 

FE  ( 1) 

10.0 

24.0 

Transverse  Pressure  Angle 

APTD(l) 

18.0 

30.0 

Helix  Angle 

AAD(l) 

10.0 

40.0 

Pinion  Base  Diameter 

DBP(l) 

8.0 

40.0 

Active  Tooth  Depth 

ATD(l) 

5.0 

40.0 
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TABLE  20 


ADS  PARAMETERS 

ADS  Parameter 

NRA 

NCOLA 

NRWK 

NRIWK 

INFO 

ISTRAT 

IOPT 

IONED 

IPRINT 

IGRAD 

NDV 

NCON 

X(l) 

X(2) 

X(3) 

X(4) 

X  ( 5) 

VLB (1) 

VLB (2) 

VLB(3) 

VLB (4) 

VLB ( 5) 

VUB(l) 

VUB ( 2) 

VUB (3) 

VUB (4) 

VUB (5) 

OBJ 


FOR  NOISE  MINIMIZATION 

_ MARGO  Value 

10 

5 

1000 

800 

0 

0 

1 

9 

2220 

0 

5 

0 

FE  (1) 

APTD(l) 

AAD(l) 

DBP(l) 

ATD(l) 

10.0 

18.0 

10.0 

8.0 

5.0 

24.0 

30.0 

40.0 

40.0 

40.0 

ACR(l)  (Aggregate  Contact 
Ratio) 
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